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ABSTRACT 


An  investigation  was  conducted  on  the  distribution  and  magnitudes 
of  boundary  shear  stresses  arising  from  subcritical  flows  through  curved 
trapezoidal  channels.     The  test  series  was  designed  to  determine  the 
effects  on  the  shear  patterns  of  variation  in  dis charge ,  bend  geometry, 
boundary  roughness ,  and  upstream  channel  alignment.    Two  flumes  of 
different  bend  radius  and  base  width  were  used,  each  consisting  of  a 
single  circular  curve  of  60°  central  angle,  with  straight  upstream  and 
downstream  tangent  sections.    Both  flumes  were  constructed  with  2  to  1 
side  slopes,  and  the  base  horizontal  along  all  transverse  sections.  The 
greatest  part  of  the  test  program  dealt  with  flows  in  smooth  channels. 
For  a  Froude  number  range  of  0.32  to  0„£5,  the  stream  geometries  varied 
as  follows?    ratio  of  width  to  depth,  7  <  w/y0  <  12|  ratio  of  width  to 
centerline  radius,  0.29  <  w/rc  <  0,80,    Two  tests  were  conducted  in  a 
rough-surfaced  channel  at  stream  geometries  corresponding  to  runs  in  the 
smooth  channel  series „    For  the  final  four  tests  described,  which  were 
conducted  in  a  smooth  channel,  sets  of  screens  were  installed  in  the 
approach  flow  in  order  to  simulate  the  disturbances  caused  by  additional 
curves  upstream  from  the  test  reach.    By  this  method,  the  shear  distri- 
bution was  studied  for  systems  of  two  curves  of  similar  sense,  and  for 
sets  of  reverse  curves.    Local  boundary  shear  stresses  were  measured 
with  surface  Pi tot  tubes  adapted  for  application  in  free  surface  flows. 
The  calibration  for  these  instruments,  originally  developed  for  air  flows 
through  smooth  pipes,  was  found  to  be  valid  for  direct  application  in 
the  smooth  channels .    A  modified  Pi tot  tube  was  developed  and  calibrated 
for  use  on  the  rough  test  surface , 

The  boundary  shear  stress  data  are  presented  as  relative  shear, 
^o/^o  (i.e.  as  local  shear  in  terms  of  the  shear  for  uniform  flow), in 
the  form  of  contour  maps  of  the  test  reach|  velocities  and  water  surface 
elevations  are  presented  by  section. 

The  boundary  shear  patterns  obtained  cannot  be  predicted  quanti- 
tatively from  the  gross  characteristics  of  the  flow.    Local  shears  were 
found  to  occur  at  intensities  of  more  than  double  the  mean  tractive 
force  computed  for  uniform  flow  3  as  might  be  expected,  the  intensities 
of  these  local  shears  increase  markedly  with  the  stream  curvature.  A 
one-dimensional  treatment  of  the  energy  dissipation  in  the  curve  and 
downstream  tangent  fails  to  indicate  either  the  peripheral  locations  or 
the  intensities  of  the  greatest  boundary  stresses,,    In  general,  the 
distributions  and  the  relative  magnitudes  of  the  local  boundary  shear 
stresses  appear  to  be  functions  of  the  stream  geometry,  for  streams  in 
subcritical  motion. 

The  patterns  of  velocity  and  shear  in  a  stream  curve  are 
influenced,  both  directly  and  indirectly,  by  the  transverse  and  longi- 
tudinal pressure  gradients  in  the  curve.    An  analysis  of  the  water  surfac 
superelevation  in  curving  flow  reveals  that,  for  conditions  of  moderate 
curvature,  the  transverse  water  slope  is  quite  insensitive  to  variations 
in  the  velocity  distribution,  and  that  it  depends  only  on  the  mean 
momentum  of  the  flow  and  the  stream  geometry. 


ACKNOWLEDGEMENTS 


This  investigation  was  conducted  at  the  Hydrodynamics  Laboratory , 
Department  of  Civil  Engineering,  of  the  Massachusetts  Institute  of 
Technology  during  the  period  March  1958  through  September  1961.    The  work 
was  sponsored  by  the  Soil  and  Fater  Conservation  Research  Division, 
Agricultural  Research  Service,  U.  S.  Department  of  Agriculture,  under 
Contract  Number  12-lU-100-2£90[Ul] ,  and  Cooperative  Research  Agreement 
Number  12-lU-100-5>227 [hi] .    The  support  for  this  work  was  established 
through  the  MIT  Division  of  Sponsored  Research  under  DSR  Project  Numbers 
5-7922  and  U-8682. 

The  experimental  study  was  outlined  in  close  liaison  with  Messrs. 
Louis  M.  Glymph,  Jr.,  Carl  R.  Miller,  and  Donald  A.  Parsons  of  the 
Agricultural  Research  Service.    Their  interest  in  the  work  contributed 
greatly  to  its  success,  and  through  periodic  re-evaluation  of  the  project 
aims,  new  areas  of  study  were  opened  which  had  not  been  envisioned  in  the 
original  proposals.    Dr.  Arthur  T.  Ippen  acted  as  project  supervisor 
throughout  the  investigation.    Dr.  Ronald  E.  Nece,  and  Messrs.  Charles  A. 
Givler  and  Georges  K.  Noutsopoulos  participated  in  the  early  phases  of 
the  work.    Messrs.  William  R.  Jobin  and  Omar  H.  Shemdin,  Research 
Assistants,  performed  much  of  the  experimental  work  and  prepared  the 
drawings  for  the  report.    Dr.  Philip  A.  Drinker  was  associated  with  the 
project  throughout,  initially  as  Research  Assistant,  and  during  the 
final  year  of  study  as  Geologist,  Soil  and  ¥ater  Conservation  Research 
Division,  Agricultural  Research  Service. 


TABLE  OF  CONTENTS 

Page 

I.      INTRODUCTION  1 

II.      FLOW  THROUGH  CURVED  CHANNELS:     A  REVIEW  2 

III.      EXPERIMENTAL  EQUIPMENT  AND  TECHNIQUES  7 

A.  The  Test  Channel  7 

1.  Description  of  the  Equipment  7 

2.  Operating  Conditions  8 

B.  Instrumentation  8 

1.  The  Measurement  of  Boundary  Shear  Stress  ;8 

a.  General  Requirements  8 

b.  Shear  Measurements  by  Surface  Pi  tot  Tubes  '.  9 

c.  Experimental  Application  11 

i.  Smooth  Surface  H 

ii.  Rough  Surface  11 

2.  Additional  Instrumentation  12 

a.  Manometry  12 

b.  Depth  Measurements  12 

c.  Velocity  Measurements  12 

IV.      THE  EXPERIMENTAL  PRO (31  AM  1$ 

A.  Scope  of  the  Study  l£ 

B.  Presentation  of  Data  16 

1.  Velocity  and  Water  Surface  Profile  Data  16 

2.  Boundary  Shear  Stress  Data  17 

V.      DISCUSSION  OF  RESULTS  20 

A.  Tests  Conducted  with  Uniform  Conditions  of  Approach: 

Flow  Through  Single  Curves  20 

1.  Smooth  Channel }  Runs  1-7  20 

a.  Entrance  Conditions  20 

b.  Flow  through  the  Curve  '20 

c.  Boundary  Shear  Stress  2k 

2.  Rough  Channel ,  Runs  8  and  9  35 

a.  Entrance  Conditions  35 

b.  Flow  through  the  Curve  35 

c.  Boundary  Shear  Stress  35 

B.  Tests  with  Simulated  Compound  Curve  Systems: 

Runs  10-13                        '  38 

1.  Entrance  Conditions  3$ 

2.  Flow  through  the  Curve  39 

3.  Boundary  Shear  Stress  hY 


-iii- 


TABLE  OF  CONTENTS  (continued) 

Page 

C.  Superelevation  in  the  Channel  Bend  H8 

D.  Energy  Dissipation  in  a  Channel  Bend  Jii- 

E.  Scour  Patterns  in  Alluvial  Streams  6U 

VI.      CONCLUSIONS  67 

A.  Flow  Through  Single  Curves  67 

B.  Compound  Curve  Systems  68 

C.  General  Observations  68 

VII.      BIBLIOGRAPHY  70 

VIII.      APPENDIX  -  Superelevation  in  Curved  Flow  7U 


-iv- 


LIST  OF  FIGURES 

Figure  Page 

1.  Zones  of  Separation  in  a  Laboratory  Channel  Bend  6 

2.  Helicoidal  Motion  in  a  Channel  Bend  6 

3.  The  Two  Trapezoidal  Test  Channels  13 

U.      Model  III  Surface  Pi tot  Tube  lh  ' 

Velocity  Distributions;    Run  No.  1  22 

6.  Velocity  Distributions:    Run  No.  1|-B  22 

7.  Velocity  Distributions?    Run  No.  6  22 

8.  Transverse  Fater  Surface  Profiles,  all  Stations, 

Run  No.  1  '  23 

9.  Transverse  Fater  Surface  Profiles,  Station  5>, 

Runs  No.  5,  6,  7  23 

10.  Peripheral  Distribution  of  Boundary  Shear  in  Straight 
Trapezoidal  Channels  26 

11.  Boundary  Shear  Contour  Map;    Run  No.  1  27 

12.  Boundary  Shear  Contour  Map;    Run  No*  2  28 

13.  Boundary  Shear  Contour  Map;  Run  No.  3  29' 
ill.  Boundary  Shear  Contour  Map;  Run  No.  I4.  30 
1$.      Boundary  Shear  Contour  Map;    Run  No.  5  31 

16.  Boundary  Shear  Contour  Map;    Run  No.  6  3.2 

17.  Boundary  Shear  Contour  Map;    Run  No.  7  33 

18.  Velocity  Distributions:    Run  No.  9  36 

19.  Transverse  Fater  Surface  Profiles:    Run  No.  9  36 

20.  Boundary  Shear  Contour  Map:    Run  No.  8  37 

21.  Boundary  Shear  Contour  Map:    Run  No.  9  h0 

22.  Velocity  and  Shear  Distributions  in  the  Approach  Flow: 

Run  No .  10  '  la 

23.  Velocity  and  Shear  Distributions  in  the  Approach  Flow: 

Run  No.  11  111 

2)4..      Velocity  and  Shear  Distributions  in  the  Approach  Flow: 

Run  No .  12  '  la 

25.  Velocity  and  Shear  Distributions  in  the  Approach  Flow: 

Run  No .  13  la 

26.  Transverse  Fater  Surface  Profiles:    Run  No.  11  1|2 

27.  Transverse  Fater  Surface  Profiles:    Run  No.  13  b2 


-v- 


LIST  OF  FIGURES  (continued) 

Figure  Page 

28.  Boundary  Shear  Contour  Maps    Run  No.  10  k3 

29.  Boundary  Shear  Contour  Map:    Run  No.  11  hh 

30.  Boundary  Shear  Contour  Map:    Run  No.  12  U5 

31.  Boundary  Shear  Contour  Map:    Run  No.  13  lj.6 

32.  Plot  of  Computed  Maximum  Superelevation  £2 

33.  Shear  Distribution  by  Sections:    Runs  No.  1  and  U-B  5? 

3U.      Average  Rate  of  Energy  Dissipation  through  the  Curve: 

Runs  No.  1  and  U-B  58 

35.  Variation  with  Curvature  of  the  Total  Mean  Shear  Stress  60 

36.  Variation  with  Curvature  of  the  Maximum  Relative  Shear  60 

37.  Variation  with  Curvature  of  Scour  Areas  62 

38.  Map  and  Sections  of  an  Alluvial  Stream  61 
A-l  Superelevation  Computed  for  a  Free  Vortex  78 
A- 2      Superelevation  Computed  for  a  Forced  Vortex  78 

LIST  OF  TABLES 

Table  Page 

I.      The  Range  of  Channel  Geometries  Covered  in  the  Test 

Program  18 

II.      Summary  of  Test  Conditions  Determined  in  the  Approach 

Flow    "  19 

III.      Maximum  Fater  Surface  Superelevations ,  Measured  and 

Computed  53 

IV.      Comparison  of  the  Average  Rates  of  Energy  Dissipation 

in  the  Test  Reach                                       "  59 


-vi- 


DEFINITIONS  AND  NOTATIONS 


A 

=    channel  cross-section  area. 

b 

cfx 

=    bottom  width  of  channel. 

=    mean  coefficient  of  friction  at  section  x  = 

Tx/  2 

d 

diameter  of  Pi tot  and  Prandtl  tubes . 

F 

=    Froude  number  =  V  /  \/gy^  . 

g 

=    gravitational  acceleration,  ft/sec2. 

Ho 

=    specific  head,  y0  +  V2/2g,  ft- lb/lb. 

AH 

=    manometer  reading,  ft.  of  water. 

k 

absolute  height  of  roughness  particle. 

ke 

=    equivalent  sand  roughness  height. 

L 

=    center line  arc  length  of  the  test  reach. 

Pt 

-    total,  or  stagnation,  pressure,  lb/ft2. 

Po 

=    static  pressure,  lb/ft2. 

P 

. =    wetted  perimeter. 

Q 

=    volumetric  rate  of  flow,  ft3/sec. 

r 

-    radial  distance  to  a  point  in  curve. 

rc 

=    radius  of  curvature  of  channel  center line. 

ro 

=    radius  of  curvature  of  outer  edge  of  water 

surface . 

R 

=    hvdraulic  radius.  A/P.  ft. 

E 

=    Reynolds  number,  IiRV/v, 

=    boundary  layer  Reynolds  number  Vx/v. 

% 

=  shear  Reynolds  number,  u^y/v. 

So 

=    channel  slope,  ft/ft. 
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DEFINITIONS  AND  NOTATIONS  (continued) 


Se  f 

energy  gradient  ft/ft. 

=    average  energy  gradient  over  the  test  reach. 

S*x  = 

=    local  rate  of  energy  dissipation  at  section  x,  ft-lb/lb/ft. 

T  = 

=    water  temperature ,  °F. 

u,v  = 

3    local  velocity,  ft/sec. 

A 

•    shear  velocity,    \/t  /p,  ft/sec. 

V  = 

=    free  stream  velocity  or  average  velocity,  Q/A,  ft/sec. 

w  = 

-    water  surface  width. 

x  = 

=    longitudinal  distance  measured  along  channel  centerline . 

y  = 

=    distance  normal  to  plane  of  channel  bottom. 

=    average  depth  =  A/w. 

y0  = 

=    total  depth  of  flow. 

z  = 

=    transverse  distance. 

a  = 

-    angular  misalignment  of  test  probe. 

T  = 

-    specific  weight  of  liquid,  lb/ft3. 

5  * 

:    boundary  layer  thickness,  ft. 

s  = 

:    kinetic  energy  correction  factor. 

9  = 

■    local  central  angle  of  curvature  from  beginning  of  bend,  degrees. 

-    dynamic  viscosity  of  fluid,  lb-sec/ft2. 

v  = 

8    kinematic  viscosity  of  fluid,  ft2/sec. 

P  s 

=    mass  density  of  fluid,  slugs/ft3. 

To  = 

■■    local  boundary  shear  stress,  lb /ft2. 

"o  = 

=    average  boundary  shear  stress  at  curve  entrance (Sta. l)lb/f t2 . 

Tx  = 

:    average  boundary  shear  stress  at  section  x. 

T  ■ 

"    average  boundary  shear  stress  over  the  full  test  reach. 

f (),$()    =    function  of. 
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I.  INTRODUCTION 


The  destruction  of  valuable  land  by  erosion  in  stream  bends  is  a 
problem  which  has  stimulated  a  considerable  amount  of  research  directed 
towards  the  control  of  meandering  and  scour  in  alluvial  streams.    As  an 
immediate  attack  on  the  problem,  primary  interest  has  been  focused  on 
techniques  of  channel  stabilization  through  the  placement  of  erosion- 
resistant  revetments  or  other  boundary  modifications  designed  to  inhibit 
local  scour.    For  lack  of  systematic  quantitative  information  on  the  modes 
of  erosion  in  curves,  the  installation  of  revetments  has  an  empirical 
basis,  founded  upon  field  studies  of  localized  stream  systems .  Laboratory 
studies  of  sediment  motion,  as  well  as  of  velocity  patterns  in  curving 
flumes,  have  contributed  greatly  to  the  general  understanding  of  stream 
hydraulics  and  morphology.    However,  owing  to  the  complexity  of  the 
motion  of  fluids  in  curved  channels,  it  has  not  been  possible  either  to 
establish  a  satisfactory  theoretical  basis  for  curving  flow,  or,  failing 
this,  to  predict  the  scour  intensities  from  measured  velocity  distributions. 

The  study  described  in  this  report  was  conducted  to  determine  the 
distributions  and  magnitudes  of  the  boundary  shear  stresses  generated  in 
curved  reaches  of  rigid  prismatic  channels.    In  this  work,  which  represents 
an  initial  attack  on  the  erosion  problem,  all  sediment  properties  were 
excluded  from  consideration,  and  only  the  effects  of  the  flow  pattern  on 
a  clearly  defined  set  of  boundary  surfaces  were  studied.    A  general 
explanation  of  the  various  related  aspects  of  sediment  mechanics  — 
erosion,  transport,  and  deposition  —  will  require  much  supplementary 
investigation  of  systems  containing  sediment  in  order  to  establish  first, 
the  role  of  sediment  in  modifying  the  dynamic  properties  of  the  trans- 
porting fluid  itself,  and  second,  the  effect  on  the  scour  pattern  of 
boundary  modifications  caused  by  the  movement  of  bed  material. 

An  earlier  report  (19)  describes  the  development  of  the  test 
equipment  and  instrumentation,  and  presents  the  results  of  the  tests 
conducted  during  the  period  covered  by  the  original  research  program. 
In  order  to  broaden  the  scope  of  the  investigation,  a  further  series  of 
tests  was  performed  in  which  the  range  of  stream  parameters  was  extended. 
This  report  comprises  the  final  report  on  these  studies  of  boundary  shears 
in  trapezoidal  channels.    In  order  to  better  coordinate  the  analyses  and 
discussions,  the  results  of  all  tests,  including  those  recorded  previously, 
are  presented.    The  experimental  equipment  and  instrumentation  have  been 
described  at  length  (19)  and,  except  for  the  modifications  in  the  test 
channel,  equipment  details  are  not  included  here. 

In  view  of  the  findings  from  the  later  experiments ,  the  original 
data  were  re-examined,  resulting,  in  several  cases,  in  modifications  of 
ideas  formed  earlier.     In  substance  however,  the  conclusions  drawn  from 
the  earlier  tests  were  not  greatly  affected  by  the  results  of  the  final 
experiments;  rather,  it  will  be  seen  that  the  increased  scope  of  the 
investigation  lends  greater  generality  to  these  conclusions a 
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II.    FLOW  THROUGH  CURVED  CHANNELS:     A  REVIEW 


As  a  basic  model  of  flow  through  channel  curves ,  consider  first 
the  flow  through  a  deep,  smooth-surfaced,  prismatic  bend.    If  frictional 
effects  are  assumed  negligible,  potential  theory  permits  description  of 
the  curved  flow  as  a  free  vortex.    Mockmore  (31),  among  others,  has 
demonstrated  that  at  some  distance  from  the  boundaries  both  the  velocity 
distribution  and  the  transverse  water  surface  profile  are  well  described 
by  an  irrotational  vortex.    For  this  case  the  transverse  velocity 
distribution  is  given  by  the  constant  product  of  local  tangential 
velocity,  u,  and  radius,  r.    By  equating  the  radial  forces  on  a  fluid 
particle  as, 


.2 


■  p  - ,  (i) 

dr  r  3 

the  water  surface  profile  can  be  computed  from 

Ay  =  ijr dr  >  (2) 


where  the  flow  in  the  curve  is  assumed  to  be  concentric.    For  the  free 
vortex  the  water  surface  profile  is  hyperbolic,  convex  upward. 

However,  examination  of  the  flow  near  the  boundaries  shows  that 
even  in  a  smooth-walled  flume,  frictional  effects  combined  with  the 
superelevation  tend  to  cause  two  important  deviations  from  potential  flow: 
helicoidal  motion  of  the  main  stream  and  separation  in  zones  of  adverse 
pressure  gradient. 

In  1876  Thompson  (hh)  gave  his  now  familiar  explanation  of 
helicoidal,  or  spiral  flow  as  being  due  to  the  difference  in  centripetal 
acceleration  between  the  fluid  near  the  bottom,  which  has  been  retarded 
by  boundary  drag,  and  the  faster  moving  layers  near  the  free  surface. 
This  acceleration,  which  results  from  the  radial  increase  in  depth  (and 
hence  in  pressure),  causes  the  low  velocity  fluid  near  the  bed  to  move 
towards  the  center  of  the  curve,  while  the  fluid  near  the  surface,  having 
excess  momentum,  moves  towards  the  outside. 

The  superelevation  of  the  water  surface  in  a  channel  curve  leads 
to  two  occurrences  of  longitudinal  increase  in  depth:    along  the  outer 
bank  of  the  curve  as  the  water  rises  entering  the  curve,  and  at  the  curve 
exit  along  the  inner  bank  where  the  flow  recovers  to  normal  depth  (see 
for  example,  Figures  5>  and  8).    For  channels  of  large  curvature,  separation 
may  develop  in  these  two  areas  and  the  water  surface  may  slope  upstream 
locally.     The  fluid  particles  near  the  boundary  have  insufficient 
kinetic  energy  to  move  against  a  pressure  gradient,  and  separation  of  the 
boundary  layer  from  the  streambank  results.    A  simple  circular  curve 
connecting  two  tangent  channels  is  likely  to  produce  larger  separation 
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zones  than  would  a  natural  stream  bend,  with  more  gradual  transitions. 
In  periods  of  high  flows  separation  in  natural  channel  bends  manifests 
itself  as  eddies  or  "bank  rollers"  and  results  in  bar  deposits. 

It  has  been  pointed  out  (26)(U6)  that  superelevation  as  such 
in  subcritical    curving  flow  is  quite  insensitive  to  radial  variation 
in  the  velocity  distribution.    Substitution  into  Equation  (2)  of  other 
assumed  velocity  distributions ,  and  integration  over  the  full  width  of 
the  stream,  give    only  minor  variations  in  the  total  elevation  differ- 
ence, even  though  the  surface  profiles  vary  in  shape.    Superelevation  is 
a  result  of  curved  flow  for  all  fluids,  and  is  essentially  independent 
of  the  frictional  aspects  and  the  velocity  distribution  of  the  flow. 
It  is  controlled  primarily  by  the  boundary  geometry  and  the  mean  momentum 
of  the  stream.    Thus,  while  superelevation  provides  the  mechanism  by 
which  separation  and  helicoidal  motion  take  place,  their  direct  cause  is 
boundary  friction  or  shear  stress,  and  in  a  completely  frictionless  flow 
they  could  not  occur. 

In  deep,  turbulent  streams,  the  separation  zone  at  the  concave 
bank  is  usually  obliterated  by  the  turbulent  infusion  of  momentum  from 
the  main  body  of  the  stream,  although  bank  rollers  are  observed  at  the 
outside  of  abrupt,  shallow  bends.    The  separation  zone  at  the  inner  bank 
is  more  persistent,  as  shown  by  the  point  bar  deposits  common  to  virtually 
all  streams  except  those  which  are  deeply  incised  and  actively  degrading. 
The  persistence  and  extent  of  this  latter  zone  are  due  both  to  the  abrupt 
rise  in  water  surface  at  this  point,  and  to  the  spiral  motion  of  the 
stream  which  tends  to  move  the  faster  moving  fluid  to  the  opposite  bank. 

Examples  of  helicoidal  motion  and  separation  in  a  laboratory 
flume  bend  are  shown  in  Figures  1  and  2.    Both  of  these  photographs  were 
taken  from  the  bend  exit,  facing  upstream,  so  that  the  flow  is  from  top 
to  bottom  of  the  photographs .    Figure  1  was  taken  after  dye  had  been 
injected  and  allowed  to  spread  to  the  full  extent  of  the  separation 
zones.    In  Figure  2  the  fainter  dark  streaks  were  caused  by  dye  crystals 
moving  towards  the  inner  bank  with  the  bottom  current  (the  dashed  arrow), 
while  the  heavier  plume  was  made  by  continuous  injection  of  dye  at  the 
surface    (the  solid  arrow). 

These  two  phenomena  have  far  reaching  effects  on  the  flow  in 
natural  streams.    The  helicoidal  pattern  is,  perhaps,  the  more  important, 
in  that,  by  shifting  the  faster  moving  fluid  to  the  outside  of  the  bend, 
it  tends  to  suppress  the  free  vortex  pattern.    Einstein  and  Harder  (8) 
show  experimentally  that,  for  a  curve  of  sufficient  length,  a  stable 
velocity  distribution  develops  in  which  the  velocity  increases  radially 
outwards.     That  is,  of  course,  the  phenomenon  observed  in  long  river 
curves  in  which  the  thalweg  remains  near  to  the  concave  bank. 

The  pattern  of  helicoidal  flows  in  the  bend  of  a  rectangular 
flume  was  studied  by  Shukry  (39)  using  a  Pi  tot  sphere  to  determine  the 
three  components  of  the  local  velocity  vector.     It  was  found  that  the 
helicoidal  (i.e.  radial)  velocities  increase  with  increasing  relative 
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curvature,  w/rc,  and  ratio  of  width  to  depth,  w/y0.    Einstein  and  Harder  (8) 
also  observed  that  increasing  the  ratio,  w/y0,  as  well  as  increasing  the 
boundary  roughness,  led  to  high  velocities  towards  the  outer  bank,  due, 
presumably,  to  a  more  pronounced  helicoidal  motion. 

The  laboratory  bends  in  which  the  irrotational  type  of  flow  is 
observed  are  too  short  to  complete  the  transfer  of  the  filament  of 
highest  velocity  to  the  outer  bank.    The  patterns  of  flow  in  sharp  bends 
—  of  both  flumes  and  rivers  —  represent  a  transitional  type  which, 
with  sufficient  length  of  curved  reach,  would  develop  into  the  friction 
controlled  flow  described  above.    The  length  of  the  transition  reach, 
however,  has  not  been  well  defined  in  terms  of  the  pertinent  parameters. 

Recently  Leopold  and  co-workers  (26)  have  reported  an  investi- 
gation on  the  frictional  resistance  of  rigid  sinuous  channels.  By 
studying  the  overall  head  loss  over  a  reach,  for  channels  of  constant 
shape  but  varied  sinuosities,  the  authors  show  that  not  only  is  the  total 
energy  dissipation  greater  in  curved  channels  than  in  straight,  but  that 
it  increases  dis continuously  above  a  certain  threshold  Froude  number,  due 
to  the  formation  of  local  standing  waves,  or  "spills",  along  the  banks. 
While  the  rate  of  increase  in  total  boundary  shear  is  not  clearly  defined, 
a  relationship  is  suggested  between  the  ratio  of  width  to  radius  and  the 
Froude  number  at  which  the  discontinuity  occurs.     In  the  laboratory  tests 
the  threshold  Froude  number  varied  from  about  0 . U  to  0.6,  while,  from  a 
survey  of  river  data,  the  authors  found  that  bankfull  flows  generally 
occur  in  streams  in  the  range  0.2  <  F  <  0.U5.    They  suggest  that  the  rare 
occurrence  of  higher  Froude  numbers  may  be  due  to  a  tendency  of  alluvial 
rivers  to  eliminate  local  areas  of  high  resistance.    Such  a  local  reduction 
in  boundary  shear  could  come  about  through  erosion  of  the  banks  at  points 
of  greatest  spill  occurrence,  where  there  is  increased  dissipation  of 
energy. 

The  relationship  between  shear  stresses  and  erosion  is  further 
complicated  in  natural  streams  by  loose  bed  material  which  permits  an 
additional  degree  of  freedom  not  possessed  by  rigid  boundary  flumes. 
This  freedom  to  shape  its  own  bed  undoubtedly  causes  a  river  to  alter  its 
flow  pattern  in  ways  which  cannot  readily  be  predicted  from  a  study  of 
rigid,  prismatic  channels  (23).    Thus,  scour  will  occur  in  regions  of 
high  shear  stress;  the  modifications  in  the  bed  configuration  caused  by 
the  scour  will  in  turn  modify  the  current  pattern  and  hence  the  shear 
distribution . 

Previous  investigations  of  energy  dissipation  in  curved  flow  have 
been  restricted  to  one-dimensional  analyses  for  both  open  and  enclosed 
conduits.    As  shown  by  Anderson  (2),  the  different  experimental  techniques, 
as  well  as  the  varied  definitions  of  a  loss  coefficient,  used  by  past 
investigators  have  led  to  considerable  confusion.     Ito  (21)  has  demonstrated, 
by  a  re-analysis  of  early  data,  that  the  loss  coefficient  of  a  smooth  bend 
can  be  predicted  empirically,  if  the  measurements  and  definitions  are  given 
a  common  basis.    With  respect  to  the  present  study,  a  one-dimensional 


approach  is  not  very  revealing  in  that  it  gives  no  indication  as  to  the 
distribution  of  the  boundary  shear  stresses  that  cause  the  loss  in  bends. 
Indeed,  as  is  shown  in  the  discussion  to  follow }  the  mean  section  loss 
gives  little  indication  of  the  magnitudes  of  the  local  shears ,  and  the 
location  and  intensity  of  erosive  attack  cannot  be  predicted  by  this 
means . 


Figure  1.  (above) 
Zones  of  Separation 
in  Laboratory  Channel 
Bend. 


Figure  2.  (left) 
Helicoidal  Flow  in 
Channel  Bend. 


III.     EXPERIMENTAL  EQUIPMENT  AND  TECHNIQUES 


A.     The  Test  Channel 

1.    Description  of  the  Equipment 

The  experiments  were  conducted  in  the  two  trapezoidal  flumes 
shown  in  Figure  3.     In  both  channels  the  general  arrangement  is  the 
same,  consisting  of  a  straight  approach  section  20  feet  in  length,  a 
single  curve  of  60-degree  central  angle,  and  a  straight  10-foot  exit 
section.    The  side  slopes  are  constant  at  2  horizontal  to  1  vertical 
and  the  channels  were  constructed  with  the  invert  horizontal  along 
all  radial  sections.    A  smooth  entrance  transition  from  the  stilling 
basin  at  the  upstream  end  of  the  approach  channel  reduces  the  tendencies 
for  local  separation.    Additional  control  of  depth  is  provided  by  an 
adjustable  sluice  gate  mounted  at  the  exit  section.    Details  of  the 
basic  design  of  the  channel  facility  are  given  in  Reference  19. 

The  larger  flume  (Figure  3a)  has  a  2U-inch  base  width,  and  the 
radius  of  the  curve  centerline  is  60  inches.    For  the  smaller  flume, 
which  was  laid  in  the  bed  of  the  larger,  the  base  width  dimension  was 
halved  and  the  stream  section  was  shifted  outward  to  follow  the  concave 
bank  of  the  original  curve,  giving  a  curve  centerline  radius  of  70  inches. 

In  a  curved  trapezoidal  stream,  the  relative  curvature  --  i.e. 
the  ratio  of  surface  width  to  bend  radius  —  varies  with  depth  and  thus 
cannot  be  given  solely  in  terms  of  the  boundary  geometry.    However,  the 
limiting  values  of  relative  curvature  are  those  set  by  the  base  width. 
Thus,  the  range  of  curvatures  which  may  be  obtained  in  the  two  flumes 
is  limited  by  w/rc  >  b/rc  =  0.U0  and  w/rc  >  b/rc  =  0.17  in  the  larger 
and  smaller  channels,  respectively.     In  each  channel,  therefore,  tests 
at  several  depths  give  a  moderate  range  in  relative  curvature,  the  lower 
limit  of  which  is  fixed  by  the  boundary  curvatures,  b/rCo 

During  the  installation  of  the  original  flume  structure,  the 
centerline  slope  was  set  at  0.0006U  =  l/l£63«    The  slope  was  not 
adjusted  throughout  the  test  program;  however,  after  construction  of 
the  smaller  channel,  it  was  found  that  the  approach  slope  had  decreased 
to  0.000^5  due  to  settling  of  the  foundations. 

The  channels  were  treated  with  a  waterproof  plastic  coating 
which  dries  to  form  a  tough  elastic  film.    For  the  series  of  smooth 
boundary  tests  it  was  found  that  this  material  provides  a  hydraulically 
smooth  surface,  giving  a  mean  Manning  coefficient  of  0.010o 

For  the  rough  boundary  tests  which  were  conducted  in  the  larger 
channel,  the  surface  was  prepared  as  follows?     the  channel  was  given  two 
coats  of  varnish|  while  the  surface  was  still  wet,  roughness  particles 
were  spread  by  hand  over  the  bottom  and  sidewalls  of  the  section.  After 
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the  material  had  dried,  the  surface  was  scraped  with  a  flat  bar  to  knock 
off  any  protruding  particles,  and  then  sealed  with  a  final  coat  of 
varnish.    The  rough  surface  was  applied  over  lf>  feet  of  the  approach 
section,  through  the  curve  and  for  8  feet  in  the  exit  reach.     The  rough- 
ness particles  are  uniform,  rectangular,  lucite  parallelepipeds,  of  the 
approximate  dimensions,  0.18  x  0.10  x  0.10  inches.    The  absolute  rough- 
ness height  of  the  completed  surface  was  taken  as  0.1  inches;  for  the 
treated  channel  the  Manning  coefficient  is  0.017. 

Modification  of  the  test  facility  to  incorporate  additional 
combinations  of  curves  would  have  entailed  great  additional  expense. 
Therefore,  the  effects  of  varied  channel  alignment  were  studied  through 
simulation  of  the  flow  disturbances  caused  by  curvature  in  the  approach 
channel.    For  this  purpose,  sets  of  screen  baffles  were  installed  in 
the  upstream  tangent.     The  different  screen  combinations  were  estab- 
lished by  trial  and  error  to  produce  the  desired  flow  patterns  at  the 
entrance  to  the  curve . 

2.    Operating  Conditions 

For  each  desired  depth  condition,  the  discharge  at  uniform  flow 
was  calculated  from  the  Manning  equation.    During  the  early  tests  in 
the  smooth  channels  the  Manning  coefficient  was  assumed  as  n  =  0.009. 
As  the  test  program  continued  the  channel  resistance  could  be  checked, 
and  it  was  found  that  the  proper  value  for  the  plastic  finish  was 
n  =  0.010.    A  similar  trial  procedure  in  the  rough  channel  gave  the 
Manning  coefficient,  n  =  0.017. 

With  the  proper  discharge  flowing  in  the  channel,  the  sluice 
gate  was  employed  to  adjust  the  depth  at  the  entrance  to  the  curve. 
This  downstream  control  was  necessary  because  of  the  backwater  effects 
of  the  channel,  outfall  which,  before  its  installation,  had  been 
observed  to  extend  through  the  curve  and  well  into  the  approach  tangent. 

Flow  uniformity  had  little  significance  in  the  tests  on  the 
simulated  curve  systems  due  to  the  large  scale  modification  of  the 
flow  pattern  produced  by  the  screens.    For  these  tests,  therefore,  the 
discharges  and  sluice  gate  settings  were  those  established  for  the 
corresponding  tests  at  uniform  approach  conditions. 


B.  Instrumentation 
1.    The  Measurement  of  Boundary  Shear  Stress 

a.     General  Requirements ;     The  basic  objective  of  this  study 
involved  the  mapping  of  boundary  shear  intensities  in  the  curved  test 
reach  at  varied  conditions  of  flow.    Thus,  the  primary  requirements  of 
the  shear  instrument  to  be  used  were  that  it  be  readily  movable  and  that 
it  have  a  rapid  response. 

Because  of  the  secondary  currents  in  the  bend  region,  it  was 
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further  necessary  that  the  instrument  be  insensitive  to  moderate  varia- 
tions in  the  direction  of  flow.    Local  velocities  close  to  the  boundary 
in  a  bend  possess  lateral  as  well  as  longitudinal  components ,  and 
consequently  the  boundary  shear  stresses  are  also  skewed  with  respect 
to  the  channel  center line.     For  moderate  skewness  however  (e.g.  less 
than  20  degrees),  the  magnitudes  of  the  local  velocities  and  boundary 
shear  stresses  may  be  taken  as  equal  to  their  downstream  components. 
Thus,  local  shear  magnitudes  may  be  adequately  determined  with  an 
instrument  aligned  parallel  to  the  downstream  direction,  provided  the 
instrument  itself  is  not  subject  to  appreciable  error  arising  from 
lateral  components  of  flow. 

A  discussion  of  methods  of  determining  local  shears  is  given 
in  Ref .  19;  each  of  the  techniques  listed  was  reviewed  in  terms  of  the 
requirements  of  this  study. 

b.     Shear  Measurement  by  Surface  Pi tot  Tubes:     Of  the  various 
means  available  for  the  measurement  of  boundary  shear  stress  the  surface 
Pitot  tube  technique  developed  by  Preston  (37)  was  found  to  be  best 
suited  to  the  purposes  of  this  study.    By  this  method,  the  shear  stress 
on  a  smooth  boundary  is  computed  from  the  dynamic  pressure  registered 
by  a  round  Pitot  tube  resting  on  the  surface.    Preston  demonstrated 
that,  for  a  tube  of  sufficiently  large  diameter,  the  effects  of  the 
viscous  sublayer  become  negligible,  and  the  mean  total  pressure  over  the 
face  of  the  tube  is  dependent  only  on  the  velocity  distribution  in  the 
turbulent  boundary  layer. 

The  velocity  distribution  in  the  turbulent  boundary  layer  over 
a  smooth  surface  can  be  expressed, 


uv,y  /t 
uw  v      '      -x-     W  p 


Preston  developed  the  functional  groupings  by  which  the  shear  stress  at 
the  wall,  t0 ,  could  be  expressed  as  a  dependent  variable,  and  by  direct 
calibration  in  pipe  flows  obtained  the  equation, 

t  d2  (P+-P  )d2 


valid  within  the  range, 


(pt"P0)d2 
<  l0g   <  6*5* 


Here  p  and  v  are  the  density  and  kinematic  viscosity  of  the  fluid,  and 
(pt~Po)  i-s  ^ne  dynamic  pressure  recorded  by  a  round  Pitot  tube  of  outer 
diameter  =  d.     This  calibration  equation  depends  on  the  velocity  distri- 
bution near  the  wall  which  was  found  to  be  expressible  as  the  power  law, 
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U  J 

=  8.6l(-~)      .  (5) 

Although  Preston  developed  Equation  (U)  for  a  set  of  Pi tot  tubes  of 
constant  ratio  of  inner  to  outer  diameters ,  Hsu  (l8)  later  showed, 
both  analytically  and  by  experiment,  that,  for  a  tube  of  given  outer 
diameter,  the  internal  diameter  has  negligible  effect  on  the  Pi tot 
tube  calibration. 

Equation         was  verified  by  both  Preston  and  Hsu  on  flat 
surfaces  under  varied  pressure  gradients.    For  the  present  open 
channel  work,  the  calibration  was  first  checked  directly  in  a  tilting 
glass  flume.    In  the  curve  of  the  trapezoidal  test  channel,  however, 
direct  calibration  was  not  possible  because  of  the  non-uniformity  of 
the  flow.    In  order  to  establish  the  validity  of  Equation  (U)  for  this 
application,  velocity  profiles  in  the  flow  immediately  adjacent  to  the 
boundary  (y  <  1  cm)  were  taken  at  various  points  in  the  curve.    By  this 
method  the  existence  of  the  velocity  distribution  given  by  Equation  (5>) 
was  established  (see  Reference  19,  Appendix  A);  therefore,  the  validity 
of  Equation  (U)  could  be  assumed,  subject  to  the  restriction  that  the 
surface  Pi tot  tube  lie  wholly  within  the  established  boundary  layer 
region  of  similarity. 

For  the  rough  boundary  tests,  a  surface  Pi tot  calibration  was 
developed,  which  had  been  suggested  in  principle  by  Preston.  In 
turbulent  flow  over  a  hydraulically  rough  surface,  the  velocity  distri- 
bution being  independent  of  Reynolds  number  may  be  expressed, 

-  =  f  (©  ,  (6) 
u„       a  k    3  v  ' 

where  k  is  the  absolute  height  of  the  boundary  protrusions.     Thus,  by 
dimensional  considerations  alone,  the  expression  analogous  to  Equation 
(U)  becomes, 

T°      =  •  (7) 


On  the  non-uniform  surfaces  encountered  in  general  practices,  direct 
application  of  Equation  (7)  would  be  questionable,  dur  to  the  uncer- 
tainty of  consistent  Pi tot  tube  placement.    Furthermore,  from  analysis 
of  Nikaradse's  rough  surface  velocity  data  (33),  it  seems  clear  that 
f,  in  Equation  (7),  should  vary  with  the  value  of  y0/k. 

If,  however,  tests  are  conducted  over  a  limited  range  of  y0/k 
on  a  surface  of  sufficiently  uniform  roughness,  then  it  should  be 
possible  to  find  a  single  form  of  Equation  (7).    Direct  calibration  of 
a  round  Pi tot  tube  was  performed  in  a  tilting  flume  on  an  artificially 
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roughened  surface ,  identical  to  that  which  was  later  used  in  the  test 
channel.    The  test  conditions  were;    d  -  0.1i32n,  k  =  0.1",  2"  <  yQ  <  6", 
for  mean  velocities  ranging  up  to  £  f ps .    The  empirical  equation  obtained , 

tq  =  0,02l(pt  -Pq)  ,  (8) 

is  valid  only  for  the  above  set  of  conditions.  Since  the  calibration  was 
performed  with  only  one  Pi tot  tube,  the  general  dimensionless  expression, 
Equation  (7),  was  not  established,  and  the  restricted  form,  Equation  (8), 
is  to  be  preferred. 

c.    Experimental  Application;     (i)    Smooth  Surface;    Various  models 
of  surface  Pi tot  tubes  (19; ,  were  employed  during  the  tests.    The  final 
version,  designated  as  the  Model.  Ill  instrument,  is  shown  in  Figure  U. 
The  static  tube  is  positioned  above  the  total  head  tube  in  order  to  mini- 
mize the  effects  of  the  total  pressure  gradient  near  the  boundary  on  the 
measured  static  pressure.    Before  final  assembly  of  the  instrument  it  was 
determined  that  the  two  tubes  exert  no  mutual  interference,  and  that  the 
pressure  in  the  channel  curve  is  hydrostatic ally  distributed  over  the 
depth  of  flow.    The  two  leads  from  the  instrument  are  connected  to  a 
manometer  to  give  the  desired  dynamic  head,  (p-t,  -  P0)/y« 

To  measure  the  boundary  shear  stress  at  a  point,  the  instrument 
is  mounted  vertically  with  the  total  head  tube  aligned  parallel  to  the 
downstream  direction  and  the  tip  resting  on  the  boundary.    The  shear 
value  is  computed  from  Equation  (U)  using  the  dynamic  head  as  read  from 
a  manometer.    Because  each  test  involved  about  200  readings,  it  was 
expedient  to  re-arrange  and  plot  Equation  (k)  to  give  a  direct  graphical 
solution  of  t0. 

In  preliminary  tests  it  was  found  that  the  instrument  is  quite 
insensitive  to  moderate  misalignment  of  the  probe  axis  with  the  local 
direction  of  velocity.    For  deflection  angles  up  to  20°,  the  error  in 
measured  shear  stress  varies  essentially  as  (l  -  cos  a)  \    at  a  =  20°  the 
error  is  6  °/6,  and  at  a  =  l£°  the  error  is  about  3  °A°     Dye  traces  in 
the  channel  curve  revealed  a  maximum  angularity  of  the  flow  of  about 
20°,  which  is  confined  however  to  a  very  thin  zone  adjacent  to  the 
boundary. 

(ii)    Rough  Surface;    For  the  measurement  of  shear  stress  on  the 
rough  test  surface,  the  Model  III  tube  was  fitted  with  the  adapter  sleeve 
and  bushing  shown  in  Figure  h  which  raised  the  diameter  of  the  impact 
tube  to  O.U32  inches.    This  was  done  in  order  to  increase  the  value  of 
the  manometer  readings  obtained,  as  well  as  to  reduce  the  effects  of 
local  disturbances  due  to  the  individual  roughness  particles.  The 
procedure  for  measuring  shear  stress  with  the  modified  instrument  is 
essentially  the  same  as  the  smooth  surface  technique,  the  value  of  t0 
being  computed  in  this  case  from  Equation  (8). 
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2.    Additional  Instrumentation 


a.  Manometry:    Two  types  of  manometers  were  used  in  this  study. 
These  instruments  have  been  described  in  the  earlier  report  (19).  One 
of  these,  a  micromanometer,  consists  of  two  piezometer  columns  with 
affixed  micrometer  point  gauges.    The  other,  an  inclined  manometer , 
being  the  faster  and  more  convenient,  was  used  in  most  of  the  work. 

b.  Depth  Measurements;    Depths  of  water  in  the  channel  were 
measured  with  a  point  gauge  mounted  on  the  instrument  carriage.  In 
certain  of  the  tests  it  was  more  convenient  to  determine  depths  from 
static  pressure  measurements  which  were  correlated  with  the  center line 
depths  as  determined  from  the  point  gauge  readings. 

c.  Velocity  Measurements :    Stream  velocities  were  measured  with 
a  £/l6-inch  Prandtl  tube.    The  instrument  was  always  aligned  parallel  to 
the  downstream  component  of  flow,  and  no  attempt  was  made  to  determine 
the  lateral  velocity  components . 
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Figure  3.  The  Trapezoidal  Channels.  The  locations  of  the  ten  test  stations 
are  indicated  by  the  dashed  section  lines. 
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(a.)  Instrument  Mounted  on  Flat  Surface. 
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(b)    Construction  details  showing  the  cylindrical  sleeve  adapter  for  the 
Model  III-A  Rough  Surface  Tube. 


Figure  h.    Model  III  Surface  Pi tot  tube. 
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IV.     THE  EXPERIMENTAL  PROGRAM 


A.    Scope  of  the  Study 

The  tests  were  conducted  in  the  two  channels  shown  in  Figures  3& 
and  3b ,  covering  the  range  of  flow  geometries  which  is  summarized  in 
Tables  I  and  II. 

With  the  exception  of  Runs  10-13,  the  flow  entering  the  curve 
was  essentially  uniform  and  symmetrical  about  the  centerline.  For 
purposes  of  discussion,  therefore,  it  is  convenient  to  consider  Runs 
1-9  as  comprising  one  test  series  covering  a  broad  range  of  boundary 
geometries  for  the  simple  configuration  of  a  curve  with  straight  tangents. 
Runs  1-?  comprise  the  complete  series  of  smooth  surface,  single  curve 
tests;  for  Runs  8  and  9,  the  channel  roughness  was  increased  and  flows 
were  studied  for  stream  geometries  similar  to  those  in  Runs  2  and  k* 

In  Runs  10-13,  the  velocity  distribution  in  the  approach  flow 
was  modified  through  installation  of  various  combinations  of  baffle 
screens.    By  superposing  asymmetrical  conditions  on  the  entrance  flow 
it  was  possible  to  simulate  the  effects  of  varied  upstream  channel 
alignments.     Tests  at  two  depths  were  conducted  to  determine  the  shear 
pattern  in  a  sequent  bend  of  a  double  series  (U-conf iguration) ,  and  a 
reverse  series  (S-conf iguration)  of  curves. 

During  the  course  of  the  test  series  various  modifications  were 
made  in  the  design  of  the  surface  Pi tot  tube.    The  classification  of 
instrument  models,  as  listed  in  Table  II,  is  given  with  a  description 
of  the  instruments  in  Reference  (19).    In  order  to  establish  that  the 
measured  distributions  of  relative  shear  were  not  influenced  by  the 
specific  model  designs,  the  shear  measurements  were  completely  re-taken 
with  different  instruments.    Comparison  of  the  shear  distributions  for 
Run  2-A  (Model  I)  with  2-D(Model  II),  and  Run  3-A  (Model  III)  with 
3-B  (Model  II),  showed  that,  while  the  absolute  magnitudes  of  the  local 
shears  varied  with  the  model  used,  there  were  no  essential  differences 
in  the  maps  of  relative  shear  given  in  terms  of  the  mean  shear  stress 
measured  at  Station  1,  t0. 

In  performing  the  runs  on  the  simulated  curve  systems  with  the 
Model  III  tube,  it  was  simply  necessary,  on  the  basis  of  the  above 
tests,  to  re-measure  the  shear  distribution  at  Station  1  and  then  to 
apply  the  relative  shear  data  as  computed  previously.     Thus  the  second 
value  of  t0  given  in  Table  II  for  Run  2-A  and  for  Run  U-A,  represents 
the  measurement  with  the  Model  III  instrument,  which  was  used  in  the 
corresponding  compound  curve  tests. 

It  will  also  be  noticed  that  while  two  tests  were  conducted  at 
a  6-inch  depth,  the  energy  gradient,  mean  velocity,  and  hence  the  shear 
stress  of  the  approach  flow  were  varied  markedly  between  the  two  runs. 


Unlike  Runs  2  and  3,  which  were  repeated  at  similar  conditions  in  order 
to  verify  the  experimental  techniques ,  Runs  k~A  and  U-B  represent 
essentially  different  tests.    The  purpose  of  these  tests,  which  are 
discussed  in  a  later  section,  was  to  study  the  effect  of  approach  flow 
n on -uniformity  on  the  shear  pattern  in  the  curve. 

For  all  of  the  tests,  boundary  shear  stresses  were  measured  at 
the  ten  test  stations  indicated  in  Figures  3^  and  3b.    Highly  non- 
uniform shear  patterns  in  several  of  the  runs  made  it  necessary  to  add 
two  intermediate  sections  between  Stations  7  and  8.    Shear  measurements 
were  taken  at  2,!  intervals  across  each  station,  covering  the  full  wetted 
perimeter  except  the  outer  U"  along  each  side. 

Complete  measurements  of  water  surface  elevations  and  local 
velocities  were  made  at  each  station  in  Runs  1-h*    As  the  test  program 
progressed  coverage  of  these  variables  was  decreased  in  order  to  limit 
the  amount  of  data  to  be  processed.    At  all  depths,  however,  sufficient 
data  were  taken  to  indicate  the  gross  characteristics  of  water  surface 
superelevation  and  velocity  distribution. 

It  should  be  mentioned  that  the  numbering  system  of  the  test 
runs  shown  in  Table  I  was  established  to  give  a  logical  sequence  to  the 
order  of  discussion,  and  that  it  does  not  follow  the  chronological 
order  of  the  experiments.    All  of  the  tests  involving  the  larger  smooth 
channel  (Figure  3a)  were  completed  before  application  of  the  rough 
surface.    Following  the  completion  of  the  two  tests  on  the  rough  surface 
the  channel  was  modified  to  give  the  flume  arrangement  shown  in  Figure 
3b,  and  the  final  three  tests  were  those  conducted  in  the  smaller  smooth 
surfaced  channel. 


B.    Presentation  of  Data 

The  volume  of  the  basic  data  compiled  during  the  course  of  this 
study  is  such  that  direct  presentation  of  the  tabulated  material  would 
make  analysis  difficult.    Therefore,  in  order  to  facilitate  the  dis- 
cussions of  the  results,  the  data  presented  herein  have  been  reduced 
to  graphical  form.     The  complete  compilation  of  tabulated  data  has  been 
printed  and  is  available  in  the  basic  data  supplements  of  this  report. 

The  data  figures  for  the  various  tests  are  presented  in  the 
appropriate  discussion  sections  in  the  following  text.    Before  pro- 
ceeding with  the  discussions,  however,  it  is  convenient  to  define  the 
method  of  presentation  of  each  group  of  data. 

1.    Velocity  and  ¥ater  Surface  Profile  Data 

The  velocity  distributions  and  transverse  surface  profiles  are 
presented  by  stream  sections  following  the  sequence  of  test  stations  in 
the  curved  reach.    Because  of  the  repetitive  nature  of  these  patterns 
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drawings  were  prepared  only  for  representative  tests. 


2.    Boundary  Shear  Stress  Data 

The  value  of  the  mean  measured  shear  stress  at  Station  1,  t0, 
was  determined  for  each  test  by  graphical  integration  of  a  plot  of  the 
local  shear  stress ,  tq,  over  the  wetted  perimeter; 


The  maps  of  boundary  shear  distribution  for  the  full  test  series  are 
presented  in  terms  of  relative  local  shear,  t0/  t0.     This  form  is 
preferable  to  reporting  absolute  values 9  since  application  of  the  results 
is  simpler,  and  the  shear  patterns  for  varying  flow  conditions  are 
readily  compared.    The  absolute  magnitudes  of  local  shears  may  be  com- 
puted from  the  shear  maps  using  the  values  of  t0  (Table  I)  obtained  from 
Equation  (9). 


(9) 


TABLE  I 

Summary  of  Channel  Geometries  Covered  in  the  Test  Program. 


Run  No. 

Nominal  Stream 
Dimensions 

Channel 
Dimensions 

Boundary 
Roughness 

Velocity  Distribution 
in  Approach  Flow 

w  w 

^o  c 

1 

12  o.6o 

— < 

2 
3 

10  0.67 
8.8  0.73 

b/r    =  0.40 
Fig.  3a 

Smooth 

k 

5 
6 

8  0.80 
10  0.29 
8  0.3^ 

rb/rc  =0.17 

v-  n  =  0.010 

Symmetric; 
Single  curve 
*"  with  straight 

7 

7  o.ko 

-j 

Fig.  3b 

J 

8 

10  0.67 

- 

Rough 

9 

8  0.80 

n  =  0.017 

10 

10  O.67 

^b/r,  =  O.fco 

11 
12 

13 

8  0.80 
10  O.67 
8  0.80 

Fig.  3a 

Smooth 

> 

n  =  0.010 

Asymmetric; 
s.  Flow  pattern 
modified  to 
simulate  up- 
stream curvature. 

1 

1 

1 
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V.     DISCUSSION  OF  RESULTS 


A.     Tests  Conducted  with  Uniform  Conditions 
of  Approach:     Flow  Through  Single  Curves 

1.     Smooth  Channel,  Runs  1-7 

a.    Entrance  Conditions;    For  the  series  of  tests  conducted  in 
the  smooth  channels  with  undisturbed  approach  conditions,  the  velocity , 
water  surface  and  shear  stress  configurations  at  the  entrance  to  the 
test  section  are  essentially  symmetrical  about  the  center line  axis 
(Figures  5-10).    At  conditions  of  greater  relative  curvature  and  depth 
a  slight  asymmetry  becomes  evident  at  Station  1  due  to  the  backwater 
effects  of  the  curve.     It  seems  clear  however  that  the  flow  pattern  in 
the  bend  cannot  be  ascribed  to  upstream  distortion  of  the  flow. 


The  existence  of  a  fully  developed  turbulent  boundary  layer  was 
desired  at  the  entrance  to  the  curve.    The  boundary  layer  thickness,  5, 
may  be  approximated  by  the  Blasius  expression,  assuming  turbulent  con- 
ditions from  the  channel  entrance,  x-0; 


x 


0.38 


(10) 


Substitution  of  the  values,  V=1.50  fps  and  v  -  10~£  ft2/sec,  gives 
5  inches  at  x  -  20  feet.     Because  of  the  relatively  small  variation 

in  Reynolds  numbers,  the  boundary  layer  thickness  varied  only  a  small 
amount  over  the  full  set  of  flow  conditions.     Thus  the  boundary  layer 
was  fully  developed  in  Runs  1-2  and  5-7 ,  and  essentially  so  for  the 
5-inch  and  6-inch  depths,  Runs  3  and  I4.. 

Because  the  low  channel  slope  made  it  impossible  to  obtain  the 
energy  gradient  accurately,  it  is  not  known  to  what  extent  fully  uni- 
form flow  was  achieved.    However,  it  is  believed  that  this  condition  is 
not  of  great  consequence  for  the  results  of  this  study.     It  will,  be 
noted  that  Runs  U-A  and  -B  represent  the  largest  variation  in   F  and 
Sej  yet  the  contours  of  t0/  to  for  these  tests  were  found  to  be  virtually 
identical ..     Thus,  while  the  magnitude  of  the  mean  shear  stress  varies 
essentially  as  the  mean  velocity  head,  the  distribution  of  relative 
shears  is  but  little  affected  by  variations  in  the  Froude  number  and 
energy  gradient  of  the  approach  flow, 

b.     Flow  through  the  Curve:     The  flow  through  the  channel  bend 
exhibits  the  expected  tendency  toward  free  vortex  motion,  and  in  the 
upstream  portion  of  the  curve  (Stations  3~5)  "the  velocity  patterns  and 
water  surface  configurations  are  in  general  agreement  with  those  found 
by  previous  investigators  (3l)(39).    For  the  full  range  of  conditions 
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tested  the  maximum  total  water  surface  superelevations  (See  Section  V-C 
below)  as  well  as  the  surface  profile  shapes  (Figures  8  and  9)  conform 
quite  closely  to  the  trend  predictable  for  free  vortices. 

The  irrotational  aspects  of  the  flow  are  most  pronounced  at 
conditions  of  greatest  depth  and  curvature.     In  Runs  1-h,  the  irrota- 
tional type  of  motion  persists  over  the  full  length  of  the  curve ,  and 
high  velocities  appear  along  the  outer  bank  only  at  the  start  of  the 
downstream  tangent.     (Figures  5  and  6).    However,  at  decreased  curvature 
(Figure  7),  the  transfer  of  higher  velocities  toward  the  outer  bank  is 
quite  evident  throughout  the  second  half  of  the  curve,  and  below  Station 
5,  the  irrotational  pattern  is  suppressed. 

Attention  is  called  here  once  more  to  the  phenomena  previously 
discussed  which  are  not  readily  apparent  from  the  plots  of  velocities  and 
surface  elevation:    separation  and  helicoidal  motion. 

Figures  1,  5,  and  6  indicate  the  two  zones  of  separation  that 
may  occur  at  the  points  of  locally  rising  water  surface,.    The  zones  are 
manifested  by  stable  areas  of  reversed  flow  along  the  sloping  banks. 
Dye  injected  into  the  bank  rollers  spreads  rapidly  throughout  the 
separated  regions,  but  is  slow  to  diffuse  into  the  main  stream.  The 
presence  of  these  zones  depends  primarily  on  the  relative  curvature  of 
the  stream.    At  conditions  of  lower  curvature  and  depth,  variations  in 
surface  elevation  are  more  gradual  and  are  greatly  reduced.    In  Runs 
5>-7  there  could  be  detected  only  a  vestigial  zone  on  the  inner  bank, 
and  there  was  no  sign  whatever  of  separation  along  the  outer  boundary. 

As  is  shown  in  Figure  2,  a  helicoidal  pattern,  due  to  the 
transverse  surface  gradient,  is  well  developed  in  the  smooth  channel. 
Since  lateral  components  of  velocities  were  not  measured,  it  is  not 
possible  to  evaluate  the  effects  of  bend  curvature  on  the  intensity  of 
the  spiral  currents.    However,  the  spiral  velocities  are  generated  by 
the  transverse  pressure  gradient  in  the  curve  and  it  is  to  be  expected 
that  they  should  be  subject  to  the  same  control  and  thus  will  vary  as 
the  velocity  head  and  relative  curvature.    The  angularity  of  the  flow, 
which  is  a  measure  of  the  spiral  velocities  relative  to  the  mean  down- 
stream velocity,  did  not  appear  to  vary  greatly  over  the  range  of  test 
conditions.     In  all  of  the  runs,  the  maximum  angular  deflection  which 
could  be  observed  from  dye  traces  was  about  15>°-20°  in  a  thin  zone  close 
to  the  channel  bottom;  at  l/2"  above  the  bottom,  however,  the  velocities 
appeared  to  have  little  angularity  with  respect  to  the  curve  center line. 
The  importance  of  the  secondary  bottom  currents,  as  related  to  the  shape 
of  the  inner  separation  zone  and  to  sediment  transport,  is  demonstrated 
in  Figure  2  by  the  infusion  of  dye  to  the  separation  zone  (lower  right- 
hand  corner  of  photograph)  from  the  bed  of  the  stream. 

The  ultimate  breakdown  of  the  free  vortex  patterns  and  the 
establishment  of  a  high  velocity  zone  near  the  outer  bank  are  due  to 
two  major  factors,  the  relative  importance  of  which  depends  on  the 
stream  curvature. 
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1.  The  helicoidal  motion  tends  to  move  the  high  velocity  fluid 
continuously  toward  the  outer  bank. 

2.  The  separation  zone  which  originates  at  the  inside  of  the 
curve  reduces  the  effective  area  of  the  section,  resulting  in  acceler- 
ation of  the  constricted  stream    and  a  deflection  of  the  flow  away  from 
the  inner  bank. 

In  the  tests  at  reduced  curvature,  Runs  £~7>  the  first  of  these 
factors  is  clearly  dominant.    The  separation  zone  remains  small  yet  the 
high  velocities  begin  to  develop  along  the  concave  bank  midway  through 
the  curve . 

At  higher  curvatures  (Runs  1-h)  however,  the  separation  zone  is 
strongly  developed  on  the  inner  bank,  and  the  second  effect  assumes  the 
greater  importance.     In  addition,  because  of  the  relatively  shorter 
length  of  flow  in  a  sharper  curve,  the  helicoidal  motion  is  not  as  fully 
developed,  and  the  momentum  transfer  by  this  means  is  greatly  reduced. 
This  dominant  action  of  the  separation  zone  was  demonstrated  in  the  tests 
of  the  simulated  reverse  curve  systems  (Runs  1.1  and  13),  which  are  dis- 
cussed in  a  later  section.    In  the  latter  tests,  which  were  conducted  at 
conditions  similar  to  those  of  Runs  2  and  U-A,  the  higher  velocities  in 
the  entrance  flow  were  diverted  toward  the  inner  bank,  preventing  the 
formation  of  the  separation  zone.     Dye  traces  showed  that  the  helicoidal 
motion  was  undiminished,  yet  for  these  conditions  the  high  velocity  zone 
did  not  develop  on  the  outer  bank. 

The  further  persistence  of  the  flow  asymmetry  in  a  long  straight 
reach  following  the  curve  could  not  be  determined  in  these  tests.  The 
channel  was  too  short  to  permit  a  return  to  normal  flow,  but,  from  the 
similarity  between  the  distributions  at  Stations  9  and  10,  of  both  velo- 
city and  shear,  it  is  evident  that  such  a  return  would  be  quite  gradual 
requiring  a  considerable  length  of  straight  channel,  below  the  curve. 
Shukry  (39)  found  that  uniform  flow  is  not  re-established  for  a  distance 
of  at  least  eight  stream  widths,  while  investigations  on  pipe  bends  have 
shown  that  the  rate  of  headloss  may  not  return  to  normal  within  a  length 
of  less  than  fifty  pipe  diameters  (2) (21).    The  distance  of  recovery 
following  a  curve  is  subject  to  great  variation  as  it  depends  on  the 
degree  of  non -uniformity  produced  by  the  curve  itself.     In  general  it 
may  be  surmised  that  this  recovery  length  will  be  of  the  order  of  the 
distance  required  to  establish  a  normal  boundary  layer. 

c.    Boundary  Shear  Stress:     The  distribution  of  shear  stresses 
in  the  approach  flow  is  in  general  agreement  with  the  results  reported 
by  Enger  (l.)  on  tests  performed  in  a  straight  channel,  of  similar  trape- 
zoidal section  lined  with  sand-gravel  mixtures.     In  Figure  1.0,  the  shear 
distributions  at  Station  1  are  shown  for  the  complete  series  of  tests  in 
the  single  curve.    For  comparison,  data  from  the  rough  surface  tests, 
Runs  8  and  9,  are  included,  as  well  as  from  several  of  Enger's  tests. 

The  fact  that  the  shears  measured  in  the  gravel-lined  channel 
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show  a  somewhat  greater  scatter  emphasizes  the  difficulties  encountered 
in  working  with  streams  through  non-homogeneous  materials.    While  the 
shears  measured  in  both  the  smooth  and  rough  channels  appear  quite 
consistent,  it  must  be  borne  in  mind  that,  in  natural  streams,  greater 
variation  is  to  be  expected  in  the  local  tractive,  or  scouring,  forces. 

The  close  physical  relationship  between  the  boundary  shear  and 
velocity  patterns  is  shown  by  comparison  of  Figures  %9  6,  and  7  with 
the  corresponding  shear  maps,  Figures  11,  lh,  and  16.     That  the  shear 
pattern  should  be  governed  by  the  velocity  distribution  is  to  be 
expected.    Beneath  a  zone  of  high  velocity,  the  gradient  of  velocity  and, 
hence,  the  local  rate  of  momentum  transfer  are  increased,  resulting  in 
higher  shear  stresses.    However,  the  ability  to  relate  boundary  shear 
stress  to  local  velocities  presupposes  a  rather  idealized  boundary  layer, 
which  may  be  described  by  a  single  velocity  distribution  function. 
While  such  a  law  is  realized  for  the  flow  immediately  adjacent  to  the 
boundary  (19),  nonrunif ormity  in  the  curving  stream  prevents  such 
general  formulation  for  the  full  depth'  of  flow,  and  the  relation  between 
the  shear  distribution  and  the  gross  velocity  patterns  can  only  be 
considered  on  a  qualitative  basis. 

Comparison  of  the  shear  maps,  Figures  11-17 f  reveals  certain 
trends  in  the  variations  in  shear  patterns  with  stream  geometry.  The 
maps  for  the  tests  at  greater  curvature  (Runs  1-U)  are  essentially 
similar,  showing  the  higher  stresses  developed  on  the  inner  bank 
through  most  of  the  curve    and  on  the  outer  bank  only  below  Station 
7.    With  decreasing  curvature  the  suppression  of  the  free  vortex  motion 
becomes  more  evident  in  the  shear  maps.    Thus,  while  high  shears  occur 
on  the  inner  bank  in  Run  7>  the  crossover  of  the  shear  pattern  is  de- 
veloped between  Stations  5>  and  6.     In  Run  6,  the  inner  bank  shears  are 
greatly  diminished  and  the  crossover  occurs  at  Station  £5  while  in  Run 
5,  the  asymmetry  of  the  shear  pattern  in  the  beginning  of  the  curve  is 
barely  perceptible. 

Closely  allied  with  the  above  trend  is  the  change  in  location  of 
the  maximum  local  shear  stress  with  variation  in  stream  depth  and 
curvature;  with  increasing  depth  and  curvature  the  location  of  the 
highest  shear  shifts  upstream  from  the  outer  bank  to  the  inner  part  of 
the  curve.    In  Runs        6,  7,  1*  and  2  the  maximum  shear  was  measured 
near  the  outer  bank  in  the  downstream  tangent.     In  Run  3,  there  are 
three  points  of  high  shear,  the  maximum  occurring  high  on  the  inner  bank 
at  Station       while  in  Run  I4.,  the  maximum  occurs  at  Station  "J,  with  high 
shears  extending  above  Station  3  along  the  bank. 

An  evaluation  was  made  of  the  average  relative  shear  over  the 
entire  test  reach  affected  by  the  curve,  from  Stations  1  through  9.  The 
variation  of  the  ratio  tl/  to  with  relative  curvature,  w/rc,  which  is 
shown  in  Figure  35>3  suggests  a  trend  of  increasing  average  shear  with 
increasing  curvature.     It  should  be  borne  in  mind,  however,  that  for  the 
range  of  curvatures  studied,  the  total  variation  observed  in  tl/  EZ  i-s 
of  the  same  order  of  magnitude  as  the  experimental  error  inherent  in  the 
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Figure  11.    The  Distribution  of  Boundary  Shear  Stresses  in  the  Curved  Reach  Shown 
by  Contours  of  t0/  tc;    Run  No.  1:    smooth  channel,  b/rc  «  0.1i2,  uniform  approach, 
y0  -  2.98",  w/yG  «  12,  w/rc  -  0.6,       -  0.0070  psf . 
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Figure  12.    The  Distribution  of  Boundary  Shear  Stresses  in  the  Curved  Reach  Shown 
by  Contours  of  x0/  tJ;    Run  Nos.  2-A  and  B:    smooth  channel,  b/rc  *  0.U2,  uniform 
approach,  w/y0  -  10,  w/rc  -  0.67,    (A):    y0  -  3-86",  t0  »  O.OO87  psf., 
(B):    y0  -  la. 00",  TZ  «  0.0088  psf. 
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Figure  13.    The  Distribution  of  Boundary  Shear  Stresses  in  the  Curved  Reach  Shown 
by  Contours  of  t0/        Run  Nos.  3-A  and  B:    smooth  chnnnel,  b/rc  «  0.U2,  uniform 
approach,  w/y0  «  8.8,  w/rc  -  0.73,    (A):    yc  -  5.08",       -  0.0101  psf., 
(B):    y0  -  $.08",       -  0.0121  psf.  _2?_ 
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Figure  1L.    The  Distribution  of  Boundary  Shear  Stresses  in  the  Curved  Reach  Shown 
by  Contours  of  t0/  to;  Run  Nos.  U-A  and  B:    smooth  channel,  b/rc  *  0.h2,  uniform 
aporoach,  w/y0  -  8,  w/rc  -  0.8,    (A):  y0  -  6.00",  t£  -  0.01h8  psf., 
(B):  y0  «  5.99",  ^  "  0.0060  psf. 
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Figure  IS .  The  Distribution  of  Boundary  Shear  Stresses  in  the  Curved  Reach  Shown 
by  Contours  of  t0/  to;  Run  No.  £:  smooth  channel,  b/rc  -  0.17,  uniform  approach, 
w/yQ  «  10,  w/rc  -  0.29,  y0  -  2.01",       -  0.003U  psf. 
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Figure  16.  The  Distribution  of  Boundary  Shear  Stresses  in  tne  Cirved  Reach  Shown 
by  Contours  of  t0/  to;  Run  No.  6t  smooth  channel,  b/rc  ■  0.17,  uniform  approach, 
w/y0  -  8,  w/rc  -  0.3U,  y0  -  3.02",       -  0.0055  psf . 
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Figure  17.  The  Distribution  of  Boundary  Shear  Stresses  in  the  Curved  Reach  Shown 
by  Contours  of  t0/  to;  No.  7:  smooth  channel,  b/rc  ■  0.17,  uniform  approach, 
w/yo  ■  7,  w/rc  «  0.U,  y0  -  3.98w,       -  0.007U  psf . 
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shear  measurement  technique.    Thus,  while  a  greater  range  in  w/rc  would 
be  necessary  in  order  to  define  the  variation  in  %\J  To,  the  trend  shown 
in  Figure  35>  is  clearly  indicative  of  the  effect  of  stream  curvature  on 
bend  resistance. 

The  tests  conducted  at  similar  ratios  of  width  to  depth  are 
indicated  in  Figure  3£.    For  Runs  2  and  5,  w/y0  =  10,  and  for  Runs  h  and 
6>  w/y0  =  8.    With  the  very  limited  amount  of  data  available  for  this 
comparison,  no  trend  could  be  observed  simply  on  the  basis  of  the  cross- 
sectional  shape  of  the  stream.    It  seems  reasonable  to  assume,  however, 
that  the  effect  of  depth  per  se  on  the  shear  pattern  should  be  less 
important  and  that  the  dominant  factor  is  the  stream  curvature. 

In  addition  to  the  comparison  of  the  mean  relative  shear  stresses 
shown  in  Figure  35>>  plots  were  prepared  showing  the  variation  with  w/rc 
of  the  maximum  local  relative  shear,  as  well  as  of  the  area  enclosed  by 
a  contour  of  given  relative  shear  stress.    While  the  mean  shear,  %\J  To, 
is  of  interest  as  an  indication  of  the  bend  resistance,  the  information 
given  in  Figures  36  and  37  is  perhaps  more  directly  applicable  to  the 
problem  of  revetment  design  in  channel  curves.     In  spite  of  considerable 
scatter  in  the  data,  there  is  an  evident  trend  of  increasing  relative 
shears  with  increasing  curvature. 

It  is  of  interest  to  note  that  in  the  earlier  report  (19),  which 
did  not  include  the  results  of  Runs  5-7 >  it  was  stated  that  the  mean 
relative  shear  in  Run  2  appeared  to  be  inexplicably  lower  than  those  of 
the  other  three  tests.    However,  from  the  results  of  the  completed  test 
series  (Figures  35>-37)j  it  is  evident  that  Run  2  conforms  quite  well,  to 
the  general  trend    and  that  the  relative  shears  in  Run  1  are  somewhat 
high.    The  possible  significance  of  the  deviations  of  Runs  1  and  2  from 
the  general  trend  is  discussed  in  the  section  below  on  Energy  Dissipation 
in  the  Test  Reach  (V-D). 

The  greater  deviation  of  Run  1  may  be  due  to  error  in  the  cali- 
bration of  the  static  pressure  probe  in  the  Model  II  surface  Pitot  tube, 
which  was  used  in  the  test.    At  the  time  of  the  first  four  tests,  Runs 
2  and  3  were  completely  verified  with  different  instruments,  but  since 
the  results  of  Run  1  appeared  reasonable  it  was  not  checked.  After 
completion  of  the  full  test  series,  the  remainder  of  which  was  con- 
ducted with  the  Model  III  instrument,  it  was  found  that  the  Model  II 
calibration  becomes  uncertain  at  low  depths,  tending  to  give  substantial 
positive  errors  in  the  computed  local  shears.    It  is  unfortunate  that 
this  calibration  error  was  discovered  after  the  alteration  of  the 
channel,  which  prevented  a  re- examination  of  Run  1. 

In  Figure  36  it  will  be  noticed  that  the  maximum  local  shear  in 
Run  1  conforms  to  the  indicated  trend  while  that  in  Run  2  appears  to  be 
too  low.    While  this  may  indeed  be  the  case,  it  should  be  recognized 
that  measurement  of  the  point  of  extreme  shear  is  statistically  uncertain 
for  the  relatively  small  sampling  of  the  total  stream  bed  areasj  which  was 
necessary  in  these  tests.     Thus,  a  more  reliable  criterion  of  the 
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variation  between  tests  is  the  mean  relative  shear ,  or,  alternatively, 
the  percentage  of  the  test  reach  affected  by  a  given  intensity  of 
relative  shears. 

From  these  tests  it  can  be  concluded  that  the  boundary  shear 
stresses  in  a  channel  curve  vary  with  the  flow  conditions  and,  with  the 
shear  stresses,  appear  to  be  functions  primarily  of  the  stream  geometry. 
For  the  range  of  conditions  tested  in  60-degree  curves,  7  <  w/y0  <  12 
and  1.75  <  ro/w  <  h,  the  maximum  local  shears,  expressed  as  t0/  t^,  are 
of  the  order  of  2. 

2.    Rough  Channel,    Runs  8  and  9 

a.  Entrance  Conditions:    As  with  Runs  1-7  discussed  previously, 
the  flow  at  the  entrance  to  the  curve  in  the  rough  channel  is  essentially 
symmetrical  about  the  channel  axis  (Figures  10,18,19).    While  there  is  no 
expression  analogous  to  Equation  (10)  for  the  estimation  of  the  boundary 
layer  over  a  rough  surface,  it  may  be  inferred  that  the  boundary  layer 
was  fully  developed  in  both  tests.     The  decrease  in  Reynolds  number  and 
the  increased  coefficient  of  friction  for  the  rough  channel  would  both 
lead  to  a  more  rapid  development  of  uniform  flow. 

b.  Flow  through  the  Curve;     The  tendency  towards  free  vortex 
motion  is  again  evident  in  the  velocity  distributions  and  water  surface 
configurations  for  flow  through  the  rough  channel.    Comparison  of  Figure 
18  with  Figure  6  shows  that  the  gross  characteristics  of  the  flow  are 
essentially  unchanged  by  the  introduction  of  boundary  roughness  to  a 
bend  of  a  given  geometry. 

For  channels  of  equal  cross  section  and  slope,  but  of  different 
roughnesses,  the  mean  shear  stress  over  the  full  perimeter  at  uniform 
flow,  y  RSe?  is  a  constant,  independent  of  the  roughness.    Any  variation 
in  the  friction  factors  for  different  channels  requires  an  inverse 
variation  in  the  mean  velocity  head  of  the  stream  in  order  to  maintain 
uniform  conditions.    While  it  was  not  possible  to  duplicate  exactly  the 
conditions  of  the  corresponding  smooth  channel  tests,  a  comparison  of 
the  values  of  %Z  and  y  RS-e  in  Table  II  for  Runs  2  and  8  and  Runs  U-B  and 
9  shows  that  the  magnitudes  of  the  mean  shears  at  each  depth  are  indeed 
independent  of  the  boundary  roughness. 

c.  Boundary  Shear  Stress:     The  general  patterns  of  the  boundary 
shear  distributions  for  these  tests  are  quite  similar  to  those  in  the 
corresponding  smooth  channel  tests.    As  a  result  of  the  predominant  free 
vortex  type  of  motion  the  stresses  in  the  early  part  of  the  curve  are 
greatest  on  the  convex  bank.     Towards  the  exit  of  the  curve  the  transfer 
of  high  momentum  fluid  to  the  outer  bank  causes  increased  shears  along 
the  outer  edge  of  the  stream. 

However,  a  comparison  of  Figure  20  with  Figure  12,  and  Figure  21 
with  Figure  lU  reveals  certain  deviations  from  the  smooth  channel  patterns. 
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Figure  18.    Velocity  Distributions 
by  Section,  Run  No.  9,  Rough  Channel, 
y0  -  6",  w/y0  -  8,  w/rc  -  0.8. 


Figure  19.  Transverse  Water 
Surface  Profiles  by  Section, 
Run  No.  9. 
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Figure  20.    The  Distribution  of  Boundary  Shear  Stresses  in  the  Curved  Reach  Shown 
by  Contours  of  t0/  Tq;  Run  No.  8:    rough  channel,  b/rc  *  0.h2,  uniform  approach, 
w/y0  -  10,  w/rc  -  0.67,  y0  -  3.9LW,  To"  -  0.0085  psf. 
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In  the  maps  for  the  rough  channel  tests,  the  initial  crossover  pattern 
is  quite  evident  at  Station  6,  suggesting  a  more  rapid  transfer  of  the 
stream  momentum  to  the  outer  bank.    In  addition,  while  the  shear  stresses 
throughout  most  of  the  curve  are  of  the  same  magnitudes  as  in  the  smooth 
channel  tests,  those  in  the  downstream  reach  are  considerably  higher. 

In  order  to  explain  this  behavior,  the  helicoidal  motion  in  the 
curve  may  again  be  considered.    Due  to  the  rough  boundary,  which  pre- 
vented the  motion  of  dye  crystals    and  also  caused  a  rapid  diffusion  of 
the  aqueous  dye,  it  was  not  possible  to  clearly  evaluate  the  secondary 
bottom  currents.    However,  the  presence  of  the  helicoidal  motion  was 
obvious,  and  some  inferences  can  be  drawn  regarding  its  effect  on  the 
flow  structure.    In  the  rough  channel,  the  mean  velocity  and,  hence,  the 
superelevation  in  the  curve    are  lower  than  for  the  same  curve  geometry 
in  the  smooth  channel.    Near  the  bed  of  the  rough  channel,  a  relatively 
greater  thickness  of  the  fluid  is  retarded  by  boundary  drag.    Thus,  a 
greater  portion  of  the  stream  will  undergo  transverse  motion  due  to  the 
radial  pressure  gradient.    The  result  of  this  will  be  to  accelerate  the 
crossover  bringing  more  high  momentum  fluid  to  the  outer  bank,  with  an 
ultimate  increase  in  the  attack  exerted  by  the  stream  along  the  outside 
bank. 

From  the  two  tests  in  the  rough  channel  it  may  be  concluded  that 
the  general  pattern  of  flow  through  a  bend  of  given  curvature  is  not 
greatly  affected  by  the  boundary  roughness.    With  increased  roughness  the 
helicoidal  motion  affects  more  of  the  stream,  resulting  in  a  greater 
concentration  of  momentum  and  higher  shear  stresses  along  the  outer.  For 
these  tests  at  relatively  high  curvature,  the  maximum  local  shear  stress 
is  of  the  order  of  t0/  to  =  2.5. 

B.     Tests  with  Simulated  Compound  Curve  Systems:    Runs  10-13 

1.    Entrance  Conditions 

In  the  tests  of  simulated  double  and  reverse-bend  patterns, 
the  approach  flow  was  deliberately  distorted  in  order  to  determine  the 
effect  of  upstream  curvature  on  the  boundary  shear  distribution..  For 
this  purpose,  systems  of  baffle  screens  were  installed  near  the  channel 
entrance  to  impose  the  desired  conditions  at  Station  1. 

The  simulation  of  the  double-bend  pattern  was  established  as 
follows:    using  as  models  the  downstream  shear  and  velocity  distribu- 
tions taken  in  the  corresponding  single  curve  tests  (Run  2,  Station  8 
for  Run  10;  Run  U-A,  Station  9  for  Run  11),  the  screen  sets  were 
adjusted  until  these  same  distributions  were  obtained  at  Station  1. 
Figures  22  and  23;  a,b,  and  e,      show  the  comparison  of  the  desired 
patterns  with  those  obtained  by  this  technique. 

The  reverse-bend  patterns  were  obtained  in  a  similar  way  by 
screens  placed  on  the  opposite  side  of  the  channel.    For  these  tests 


the  downstream  distributions  of  velocity  and  shear  from  Runs  2  and  I4-A 
were  reversed j  and  the  inverted  patterns  used  as  models.    The  desired  and 
obtained  distributions  are  shown  for  the  reverse-bend  tests  in  Figures  2k 
and  2£; a j  b,  and  e . 

In  each  case  the  major  portion  of  the  screen  system  was  located  as 
far  upstream  as  possible.    This  produced  a  more  gradual  shift  of  the 
velocity  distribution  near  the  entrance  to  the  curve,  without  creating 
major  surface  disturbances.     The  screen  arrangement  for  each  test  is 
shown  schematically  in  the  inset  sketch  on  the  appropriate  shear  map 
(Figures  28-31). 

The  baffle  screens  caused  non-uniform  conditions  which  prevented 
the  determination  of  approach  flow  gradient  and  which  also  increased  the 
average  shear  stress  at  Station  1  to  a  value  greater  than  that  measured 
at  conditions  of  uniform  approach.    In  the  calculation  of  t0/  tJ^,  there- 
fore, the  average  shear,  To,  was  used  as  determined  in  the  single  curve 
series  for  Runs  2  and  U-A. 

The  modifying  screen  systems  produced  velocity  and  shear  distri- 
butions corresponding  to  the  exit  conditions  of  a  curve  upstream  of  the 
channel  bend.    It  was  unnecessary  to  reproduce  a  helicoidal  motion,  since 
below  the  curve  exit  the  superelevation  and,  hence,  the  secondary  current, 
abates  rapidly.    Since  the  degree  of  similitude   that   can  be  attained 
with  simple  screen  arrangements  is  limited,  the  shear  patterns  obtained 
from  the  compound  curve  series  should  be  regarded  as  being  somewhat 
qualitative . 

2.    Flow  through  the  Curve 

The  overall  pattern  of  flow  through  the  curve  was  not  greatly 
affected  by  the  velocity  distributions  which  were  imposed  at  Station  1. 
The  helicoidal  pattern  could  again  be  detected  by  dye  traces,  although 
the  increased  turbulence  of  the  flow  caused  a  more  rapid  diffusion  of 
the  dye.    As  might  be  expected,  the  separation  zone  on  the  same  side  of 
the  flume  as  the  higher  velocities  was  reduced  while  that  on  the  opposite 
bank  was  increased  in  extent.    In  the  reverse-bend  tests,  the  roller  on 
the  convex  bank  was  barely  evident,  andj  in  the  double-bend  tests,  reduced 
velocities  occurred  on  the  concave  bank  without  reversed  flow. 

The  shapes  of  the  water  surface  profiles  (Figures  26,27)  for  the 
double  and  reverse  curves  varied  somewhat  from  those  measured  in  the 
single  curve  tests.     The  profile  at  Station  £  for  the  double  curve  test 
shows  the  effect  of  the  higher  velocities  along  the  outer  bank,  with  a 
broad  stagnant  zone  on  the  inner  part  of  the  bend.    In  the  reverse  bend 
tests,  the  separation  zone  on  the  inner  bank  is  reduced,  and  the  profiles 
are  accordingly  steeper  near  the  inside  edge  throughout  the  curve. 

In  spite  of  the  gross  variations  in  velocity  distributions  imposed 
on  the  flow,  the  maximum  superelevations  measured  in  these  tests  are 
essentially  the  same  as  those  for  the  corresponding  single  curve  tests. 
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FLOW 


Figure  21.    The  Distribution  of  Boundary  Shear  Stresses  in  the  Curved  Reach  Shown 
by  Contours  of  To/  To;  Run  No.  9:    rough  channel,  b/rc  *  0.L2,  uniform  approach, 
w/y0  -  8,  w/rc  -  0.8,  y0  -  6.05",       -  0.012  psf . 
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UNIFORM    FLOW      4"   DEPTH    LOOKING  DOWNSTREAM 

DEFINITION  1  SKETCH 

INSIDE  OF  CURVE^\  VELOCITIES  inI  FEET  PER  SECONO      ^•''OUTSIDE  OF  CURVE 


Figure  22.    Run  No.  10:  Double 
Curve,  y0  «  3.86M,  To"- 0.0081  psf, 


INSIDE  OF  CURVE 


UNIFORM    FLOW        6"  DEPTH    LOOKING  DOWNSTREAM 
DEFINITION     i  SKETCH 

VELOCITIES  IN  FEET  PER  SECONO 


OUTSIDE  OF  CURVE 


24 


STA  I  -  DOUBLE  CURVE 
SIMULATION  OF  STA  9  -ABOVE 


RELATIVE    SHEAR  COMPARISON 


r —  STA.  I  -  DOUBLE  CURVE 
-JL"    '  \  SIMULATION  OF  STA.  9-  SINGLE 
~IO^T^  CURVE 


A.  9  -  SINGLE  CURVE 


Figure  23.    Run  No._ll:  Double 
Curve,  y0  «  6.00",  t0  -  0.0123  psf. 


2    UNIFORM    FLOW   -         4    ]  DEPTH  LOOKING  DOWNSTREAM 

"*X~]i        DEFI  N  ITION  |  SKETCH 

INSIDE  OF  CURVE  VELOCITIES   IN  FEET  PER  SECOND   OUTSIDE  OF  CURVE 


24^ 

I 


STA. 8  SINGLE  CURVE 

IMAGE  REVERSED 


STA.  I    REVERSE  CURVE 

SIMULATION  OF  STATION  8-  ABOVE 


STA.  8  SINGLE  CURVE 


STA,  8  REVERSE  CURVE 


INSIDE  OF  CURVE 


OUTSIDE  OF  CURVE 


RELATIVE    SHEAR  COMPARISON 


STATION  8- SINGLE  CURVE- 
IMAGE  REVERSED 


RELATIVE  SHEAR  COMPARISON 


0.5—    STA.  I  REVERSEjCURVE 

SIMULATION  OF  STA.  9 -SI 
CURVE 


 7  _ 


0.5   ■  


-STA  9  SINGLE  CURVE 
IMAGE  REVERSED 


Figure  2U.    Run  No.  12:  Reverse 
Curve,  y0  -  3.86",  tJJ  -  0.0081  psf. 


Figure  25.    Run  No.  13:  Reverse 
Curve,  yQ  «  6.00",         0.0123  psf 


Figures  22,23,21,25.    Shear  and  Velocity  Patterns  Obtained  for  the  Simulated 
Compound  Curve  Systems.    The  four  tests  were  performed  in  the  smooth  channel, 
b/rc  «  0.h2. 
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The  superelevations  measured  and  computed  for  Runs  2,  10,  12  and  Runs 
1±-A,  11,  13  may  be  compared  in  Table  Ilia,  b. 

Because  the  velocity  distribution  in  the  single  curve  tests  was 
shown  to  be  related  to  the  mapped  shear  pattern,  the  detailed  velocity 
profiles  are  not  given  for  the  compound  curve  tests.    Parts  b  and  d  of 
Figures  22-25>  show  the  velocity  profiles  taken  at  Stations  1  and  8  in 
the  four  tests.    Further  discussion  of  the  velocity  distribution  is 
given  in  relation  to  the  shear  stress  patterns. 

3.    Boundary  Shear  Stress 

Before  proceeding  with  a  discussion  of  the  shear  stress  distri- 
butions determined  in  these  tests,  it  is  necessary  to  explain  the  shear 
maps  shown  in  Figures  28-31.     Consider  the  map  of  the  U-inch  double-bend 
test,  Figure  28;  the  right-hand  half  of  the  figure  represents  that 
portion  of  the  channel  which  was  actually  mapped  for  the  entrance  con- 
ditions of  velocity  and  shear  shown  in  b  and  e  of  Figure  22.     In  the  left 
half  of  Figure  28,  the  shear  data  for  the  single  curve  It-inch  test, 
Figure  22,  have  been  reproduced.    The  two  maps  were  joined  by  sketching 
the  contours  together  to  provide  the  shear  distribution  throughout  the 
double  bend  shown  in  Figure  28. 

For  the  U-inch  reverse-bend  test,  Figure  30,  the  right-hand 
section  was  mapped  for  the  entrance  conditions  shown  in  b  and  e  of 
Figure  2li;  the  pattern  from  Figure  12  was  reproduced  in  reverse  to 
show  the  upstream  curve.     The  maps  of  the  6-inch  double  and  reverse- 
bend  tests  were  prepared  in  exactly  the  same  way,  with  the  upstream 
curve  pattern  being  obtained  from  Figure  lU.    The  boundary  shears  were 
mapped  for  the  entrance  conditions  shown  in  Sections  b  and  e  of  Figures 
23  and  2$. 

Examination  of  Figures  28-31  reveals  the  striking  similarity 
between  the  shear  patterns  of  the  U-inch  and  6-inch  tests  for  both 
the  double  and  reverse-bend  series. 

In  the  double-bend  tests,  the  deflection  of  the  high  velocities 
to  the  outside  of  the  sequent  bend  prevents  any  degree  of  attack  on  the 
inner  bank.    In  the  downstream  reach,  comparison  of  the  contours  of 
To/  To  =         shows  that  the  shears  are  somewhat  lower  in  Run  10  than 
in  Run  11,  although  the  difference  is  slight.     The  diminution  of  the 
outer  bank  shears  within  the  sequent  curve  indicates  a  tendency  for  the 
flow  to  revert  to  the  irrotational  pattern  of  the  single  curve  tests; 
however,  the  momentum  stays  high  on  the  outside  and  leads  to  sharply 
rising  shear  stresses  on  the  outer  bank,  below  the  curve. 

Figures  30  and    31  reveal  that,  in  the  reverse-bend  tests,  high 
shear  —  higher,  in  fact,  than  any  measured  in  previous  tests  —  are 
confined  to  the  convex  bank.     This  high  shear  intensity  occurs  high  up 
on  the  bank,  throughout  the  curve,  and  the  crossover  pattern  seen  in 
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the  cerlier  tests  never  develops.    The  deflection  of  the  high  velocities 
in  the  approach  section  to  the  inner  bank  reinforces  the  natural  tendency 
towards  free  vortex  flow.    The  momentum  of  the  fluid  along  the  inner  bank 
is  sufficient  to  prevent  separation.    In  the  absence  of  the  high  velocity 
shift  to  the  outside ,  the  high  momentum  diffuses  rapidly  through  the 
stream,  and  the  shears  return  to  normal,  with  a  symmetrical  distribution 
resulting  again  in  a  short  distance. 

The  great  similarity  of  the  two  shear  patterns,  in  the  double  and 
reverse-bend  tests,  might  best  be  explained  by  the  turbulent  structure 
of  the  flow  in  the  various  tests.    For  the  single  bend  series,  the  flow 
enters  the . channel  through  a  smooth  transition  from  a  reservoir  in  which 
there  is  essentially  no  turbulence.     Thus,  at  Station  1,  all  of  the  stream 
turbulence  is  that  which  is  generated  in  the  boundary  layer.    With  the 
introduction  of  upstream  screens,  however,  free  stream  decaying  turbu- 
lence is  introduced  into  the  flow  at  the  point  of  disturbance.  Since 
the  variation  in  Reynolds  number  between  the  U-inch  and  6-inch  screen 
conditions  is  small,  it  is  to  be  expected  that  the  turbulence  intensities 
should  be  comparable.    If  the  free  stream  turbulence  were  of  sufficient 
intensity  to  overshadow  the  normal  boundary  effects,  then  the  relative 
shear  patterns  would  be  quite  similar,  as  indeed  was  the  case. 

In  each  of  Figures  22-2£,  it  will  be  noticed  that  sections  c 
and  d  are  more  closely  related  to  each  other  than  they  are  to  section  a. 
The  implication  to  be  drawn  from  this  is  that,  as  might  be  expected, 
the  effects  of  flow  disturbances  are  not  compounded  indefinitely  through 
curves . 

Within  the  limitations  of  the  simulative  approach  used  in  these 
tests,  certain  conclusions  can  be  drawn.    In  the  sequent  part  of  a  double- 
bend,  the  shear  stresses  are  highest  along  the  outer  bank  and  do  not 
vary  greatly  over  those  in  a  single  curve.    In  a  series  of  closely  spaced 
reverse  curves,  the  shears  along  the  inner  bank  may  be  greatly  increased. 
It  is  probably  safe  to  infer  that  a  tight  S-pattern,  such  as  that  in 
Runs  12  and  13,  would  not  normally  be  found  in  natural  streams;  the  high 
shears  on  the  sequent  convex  bank  indicate  that  this  configuration  would 
tend  to  straighten  itself  rapidly. 


C.     Superelevation  in  the  Channel  Bend 


As  shown  in  the  earlier  discussions,  the  major  non-ideal  aspects 
of  the  flow  through  a  curve  —  i.e.  helicoidal  motion  and  separation  — 
are  intimately  related  to  the  varied  pressure  gradients  induced  by  the 
water  surface  configuration.     It  is  of  interest,  therefore,  to  compare 
the  measured  superelevations  with  those  predicted  on  the  basis  of 
various  assumed  conditions  of  flow. 
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As  a  first  approximation  the  superelevations  may  be  computed 
from  one^dimensional  analysis.    Assuming  concentric  flow  at  a  constant 
uniform  velocity,  such  that  the  momentum  of  the  stream  acts  at  the  mean 
(center line)  radius  of  the  curve ,  the  superelevation  is  given  by, 


Ay  = 


J      gr    r      2g  r 


(11a) 


r . 
1 


where  the  subscripts ,  i  and  o,  refer  to  the  inner  and  outer  radii  of 
the  bend j  respectively. 

A  more  accurate  estimate  of  superelevation  should  be  obtained 
by  introduction  of  the  appropriate  velocity  distribution  function , 
u  =  u(r),  and  integration  of  the  resulting  surface  gradient  over  the 
width  of  the  stream.    As  is  shown  in  the  Appendix,  Section  VIII,  exact 
solutions  for  Ay  are  possible  corresponding  to  the  various  forms  assumed 
for  the  function,  u(r).    However,  because  of  the  change  in  mean  depth 
which  must  occur  within  the  curve  (see  Appendix)  these  solutions  are 
cumbersome  and  not  practical  for  direct  application.     In  the  channel 
tests,  which  were  conducted  for  a  moderate  range  of  curvatures,  it  was 
found  that  the  variation  in  measured  mean  depth  was  slight  throughout 
the  test  reach.     The  assumption  of  constant  mean  depth  greatly  simpli- 
fies the  expressions  for  superelevation,  and  in  the  following  therefore, 
it  is  assumed  that  the  mean  depth  and  hence  mean  velocity  are  constant 
over  the  length  of  the  test  reach.    The  comparison  of  the  resulting 
simplified  equations  for  Ay  to  the  general  expressions  is  given  in  the 
Appendix  discussion. 

In  these  developments  a  rectangular  channel  has  also  been 
assumed  in  order  to  permit  a  simple  definition  of  bend  geometry. 
However,  for  a  trapezoidal  stream  of  equal  width  and  mean  velocity,  the 
superelevation,  Ay,  should  be  somewhat  smaller  due  to  the  increase  in 
effective  mean  radius  as  the  flow  moves  outward  on  the  sloping  banks . 
In  addition,  for  a  real  fluid,  frictional  retardation  of  the  flow  near 
the  banks  will  further  tend  to  reduce  the  total  transverse  water  slope. 


Since  the  gross  features  of  the  flow  through  the  single  curve 
indicate  an  essentially  irrotational  behavior,  a  better  approximation  of 
the  superelevation  should  be  attainable  by  treating  the  flow  as  a  free 
vortex.    For  this  condition  ur  =  constant  and  the  total  head  is  uniform 
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and  constant  at  all  sections.    By  again  equating  the  radial  pressure 
gradient  to  the  centripetal  acceleration,  and  integrating  over  the  width 
of  the  stream,  the  superelevation,  in  terms  of  a  proportionality  constant, 
is  given  by, 

K  2 

Ay  =  4-  (-=5-  -  -nr)  ;     K  =  ur  . 

J       2g     r.^    r  2     5  1 

10 


Then,  for  uniform  specific  head,  and  by  the  assumption  of  constant  mean 
depth,  the  cross-sectional  area  is  given  by, 


r  r  ,,  2 

o  /vO  K 


A  =  yQw  =  y(rQ-r.)  =    J     y(r)dr  =    J      (Hq  - -gJ^dr  . 


r .  r . 

l  l 


This  reduces  to, 


y  =  y    =  H  -75±r —  =  H    -  tt  . 

J      ^o       o    2gr  r.        o  2g 

oi  to 

Solving  for  K]_  and  substituting  into  the  original  expression  for  Ay, 
the  superelevation  in  the  free  vortex  becomes, 


Ay  =  |!  (  1_)  (llb) 

2g  r     >_      w^  ' 


c  1- 


where 


and 


r    -  r.  =  w 
o  1 


r    +  r.  =  2r  . 
oi  c 


In  natural  streams  the  flow  in  curves  will  show  various  degrees  of 
asymmetry,  depending  on  the  entrance  conditions,  the  bend  curvature  and 
the  arc  length  of  the  curve.    In  long  stream  bends,  as  simulated  in  the 
channel  tests  by  Runs  10  and  11,  the  patterns  of  flow  ultimately  estab- 
lished within  the  curve  shows  a  radial  increase  in  velocity.     It  is  also- 
of  interest,  therefore,  to  consider  the  superelevation  produced  in  a  curve 
with  the  highest  velocities  occurring  near  the  outer  bank.    Such  a  con- 
dition may  be  approximated  by  assuming  the  stream  to  undergo  motion  as  a 
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forced  vortex,  for  which  the  velocity  increases  linearly  in  the  radial 
direction.    For  this  case, 


Av  =    2(     2_r  2)         K~  =  ■—  . 
J       2gv  o        1      '        2  r 


Again  assuming  constant  mean  depth  and  constant  average  specific  energy 
(although  the  local  specific  hear,  H0  =  y  +  u2/2g,  is  no  longer  constant), 
the  coefficient  K2  may  be  solved  from  the  expression  for  stream  area,  and 
the  superelevation  in  the  forced  vortex  simplified  to, 


Ay  = 


jr  2w/  

2§  Vl  + 


1 


12r  2 
c 


(11c) 


Figure  32  shows  a  comparison  of  the  maximum  superelevations 
computed  from  Equations  (ll).     It  is  of  interest  to  note  that,  for  the 
bends  of  moderate  curvature  normally  found  in  natural  streams, 
w/rc  ^O.^-O.?;  the  superelevation  is  virtually  independent  of  the 
radial  distribution  of  velocity.     Thus,  for  a  bend  of  given  geometry 
the  shape  of  the  transverse  surface  profile  will  vary  with  the  radial 
distribution  of  momentum;  the  total  superelevation  is  a  function  only 
of  the  mean  momentum  of  the  stream. 

For  conditions  approaching  maximum  curvature,  w/rc  =2,  (i.e. 
V±  =  0)9  the  computed  superelevation  for  a  free  vortex  becomes  infinite. 
This  is,  of  course,  in  accordance  with  the  requirement  of  infinite 
velocity  at  the  core  of  an  irrotational  vortex.     In  the  motion  of  a 
real  fluid,  separation  occurs  under  these  conditions,  imposing  a  finite 
limit  on  the  actual  magnitude  of  the  superelevation. 

The  values  of  the  superelevation  computed  for  the  complete  test 
series  from  Equations  (ll)  are  listed  in  Table  III  for  comparison  with 
the  measured  values  of  Ay.    Table  Ilia  shows  the  results  for  the  full 
width  of  the  trapezoidal  section,  while  in  Table  Illb  only  the  central 
rectangular  portion  of  the  channel,  is  treated. 

Considering  the  practical  limitations  inherent  in  the  exact 
measurement  of  water  surface  elevations,  the  agreement  between  the 
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measured  and  computed  results  is  good  for  the  full  trapezoid ,  as  well,  as 
for  the  rectangular  section.     In  the  earlier  report  on  this  study  (19), 
it  was  stated  that  the  values  measured  over  the  full  channel  width  were 
too  low.    Those  figures  were  based  on  the  measurements  made  closest  to 
the  edges  of  the  stream.     In  the  present  compilation,  shown  in  Table 
Ilia,  the  measured  profiles  were  extrapolated  to  the  stream  boundaries, 
and  the  superelevations  so  estimated  used  as  the  values  for  the  full 
stream  width.     There  is  no  real  evidence  of  the  expected  reduction  in 
Ay,  due  to  the  effects  of  the  sloping  sides,  noted  above. 

Run  1  (Table  Ilia)  shows  the  largest  deviation  from  the  predicted 
superelevation.     It  is  likely  that  in  this  case  wave  disturbances 
emanating  from  the  boundary  influenced  the  measured  value  of  Ay.  The 
possibility  of  such  error  exists  for  all  measurements  near  the  edges  of 
the  stream;  the  size  and  magnitude  of  the  error  will  depend  on  the 
distribution  of  the  disturbances  with  respect  to  location  of  the 
measurements . 

For  conditions  of  least  curvature  the  measured  superelevation 
conforms  closely  to  all  of  the  predicted  values.    At  increased  curvatures, 
in  the  single  curve  tests,  the  full-width  superelevation  predicted  for  a 
free  vortex  is  consistently  too  high,  while  the  one-dimensional  and  forced 
vortex  equations  provide  reasonably  good  estimates.     The  flow  adjacent 
to  the  inner  bank  is  retarded  by  boundary  drag  and  here  the  irrotational 
pattern  breaks  down.     It  will  be  noted  that  for  the  single  curve  tests, 
the  superelevation  within  the  rectangular  section  is,  in  general,  best 
represented  by  the  free  vortex  formulation. 

.In  the  compound  curve  tests,  Runs  10-13,  "the  superelevations  show, 
as  expected,  little  effect  of  variations  in  the  velocity  distribution. 
It  appears  that  such  differences  as  do  occur  conform  generally  to  the 
predicted  trend.     The  double  curves  (Runs  1.0  and  11.)  with  high  velocities 
on  the  outer  bank  resemble  the  forced  vortex.     In  the  reverse  curves,  the 
higher  velocities  on  the  inner  bank  reinforce  the  natural  tendency  towards 
irrotational  motion,  and  somewhat  greater  superelevation  results. 

It  may  be  concluded  that  for  streams  of  moderate  curvature  the 
water  surface  superelevation  is  quite  insensitive  to  variation  in  the 
radial  distribution  of  velocities.     In  practice  an  adequate  estimate  may 
be  obtained  through  the  simple  one-dimensional  treatment  of  the  flow, 
Equation  (lla). 


D.    Energy  Dissipation  in  a  Channel,  Bend 

In  order  to  relate  the  present  work  to  previous  studies  of  energy 
losses  in  bends,  the  measured  shear  distributions  may  be  applied  to  a 
one-dimensional,  analysis  of  the  energy  dissipation  in  the  curved  reach  of 
the  open  channel.     In  the  following  discussion  only  the  smooth  channel 
tests  for  a  single  curve,  Runs  1-7,  are  considered. 
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As  noted  earlier,  the  short  exit  tangent  made  it  impossible  to 
study  the  return  of  the  flow  below  the  curve  to  uniform  conditions.  For 
this  reason  the  results  presented  herein  will  not  bear  direct  quantitative 
comparison  to  tests  performed  in  channels  of  adequate  length.     However , 
the  detailed  information  on  the  shear  stresses  in  these  tests  contributes 
to  the  understanding  of  the  modes  of  local  dissipation  as  related  to  the 
overall  losses  in  a  bend. 

The  energy  of  a  flowing  stream  is  dissipated  through  the  resis- 
tance or  drag  upon  the  solid  boundaries.     If  a  control  volume  is  assumed 
of  cross  section  Ax, of  perimeter  Px^ and  of  length  Ax,  the  power  dissipated 
becomes : 


ART  =  t    P    -~  Ax  -  t    P    V    Ax  9 

X     X  ^*     ^  " 

wherein  tx  is  the  mean  shear  stress  over  the  perimeter  Px. 


X 


The  energy  dissipated  per  pound  of  fluid  flowing  per  unit  length  of 
channel  becomes; 


AH 
Ax 


=  S, 


x 


T  P  V 
XXX 

v  A  V 
1    x  x 


yR 

1  X 


(13a) 


For  uniform  or  normal  flow  in  straight  channels  the  quantities  involved 
become  independent  of  distance.    For  flow  through  curves  and  other 
channel  transitions,  variations  in  these  quantities  must  be  expected  with 
distance  in  accordance  with  the  non-uniform  character  of  the  flow. 


Since  the  dimensions  of  the  test-channel  for  the  experiments  on 
channel  curves  did  not  permit  the  establishment  of  completely  uniform 
flow,  the  conditions  at  Station  1  were  adopted  for  reference  purposes. 
In  order  to  avoid  confusion  with  hydraulic  properties  for  uniform  flow, 
which  are  usually  denoted  by  the  subscript,  o,  the  coordinate,  x,  is 
here  re-defined  such  that  at  Station  1,  tx  -  ^1  =  ^"o» 

Equation  (13a)  may  be  transformed  by  introducing  a  mean  friction 
coefficient, 
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Tx  =  C?xpJt  5 


AH  Tx  ?x* 

Sex  =  to  x  =  "^f"  =  °fx  2gT  '  (l3b) 

1    X  &  X 


In  these  tests,  it  was  found  that  the  maximum  increase  in  tx  is  of  the 
order  of  UO  %  while  Ax,  Px,  and  Vx  vary  only  about  2-k°fi.  From 
Equations  (13)  therefore,  neglecting  variations  in  flow  geometry,  the 
local  average  rate  of  energy  dissipation  is  proportional  to  tx.  However, 
AH/Ax|    is  only  approximately  proportional  to  Cf    within  the  variations 

of  Vx2/2g. 

In  Figure  3h  graphs  are  shown  for  two  of  the  runs  in  which  the 
rate  of  energy  dissipation  relative  to  that  at  Station  1  is  plotted  as, 

t  Se 

X  _     CX  X 

where  x  =  center line  distance  through  the  test  reach  and  rc  =  center line 
radius  of  the  bend.    The  values  of  average  relative  shear  at  each  station, 
tx/t^,  were  obtained  by  graphical  integration  of  the  section  plots  of 
shear  vs.  perimeter  shown  in  Figure  33; 

£  =  ^  =  el  (-£-)d?]x  .  (HO 

It  J  t 

o  o 

In  the  same  figure,  the  variation  in  Cf  is  plotted  as  Cf  /Cf1 . 

X  JL 

The  foregoing  relationships  may  be  applied  to  give  an  approxi- 
mation of  the  additional  head  loss  for  flow  through  the  curve.  By 
integration  of  the  tx/t^  curves,  the  average  relative  shear  stress, 
tl/t]_,  may  be  determined  over  the  test  reach  of  length,  L.  These 
computations  were  performed  for  the  test  series,  Runs  1-7 •    The  average 
relative  shears  so  obtained  are  presented  in  Table  IV  and  Figure  35>j 
which  show  the  variation  with  relative  curvature.    From  the  proportion- 
ality shown  in  Equations  (13)  between  tx  and  Sex, 

SeT  ^T 

sr-=  =        .  (15) 
Se-L  ^ 
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Figure  33.    Shear  Distribution  by  Section  Showing  Local  Values  of  t0/  t0. 
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Figure  3h.    Longitudinal  Variation  in  the  Average  Rate  of  Energy  Dissipation 
Through  the  Curve,  shown  for  Runs  1  and  U-B. 
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Figure  38.    Sketch  map  and  lateral  sections  of  an  alluvial  stream, 
showing  the  course  of  the  channel  of  maximum  scour. 
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Since  Se-^  is  the  average  energy  gradient  over  the  test  reach,  the  dif- 
ference, (Se,  -  3e<i) }  represents  the  average  increase  in  head  loss  per 
unit  length  over  the  normal  rate  of  head  loss  assumed  here  as  existing 
at  Station  1.    The  relative  average  excess  head  losses  are  shown  for 
these  tests  in  Table  IV.    Due  to  variation  in  the  location  of  Station 
10  in  several  of  the  tests,  these  computations  are  based  on  a  test 
reach  extending  from  Stations  1-9. 

It  is  appreciated  that  this  method  as  applied  here  can  give  only 
a  very  approximate  value  for  the  rate  of  energy  loss  caused  by  the  curve, 
since  the  normal  friction  slope  is  probably  higher  than  the  one  deter- 
mined for  Station  1.    Also,  the  straight  length  downstream  of  the  curve 
was  not  adequate  to  follow  the  return  of  the  shear  stresses  back  to 
normal  values.    However,  other  methods  of  determining  the  net  head  loss 
(e.g.  the  analysis  of  variation  in  specific  head)  are  subject  to  the  same 
difficulty  of  inadequate  tangent  lengths  of  channel  and  involve  the 
measurement  of  surface  gradients.    With  the  flow  conditions  employed, 
these  measurements  could  not  be  made  accurate  enough  to  give  consistent 
results,  and  hence  this  approach  of  evaluating  the  head  loss  due  to  the 
curve  was  abandoned  in  favor  of  the  direct  application  of  the  boundary 
shear  stresses. 

The  effect  of  a  curve  is  similar  to  that  produced  by  any  resis- 
tance which  might  exist  in  a  straight  channel,  such  as  a  sill,  bridge 
pier,  or  a  reach  of  increased  roughness.    Such  a  resistance  to  the  flow 
will  cause  decreased  velocities  and  increased  depths,  resulting  in  an 
upstream  storage  of  potential  energy.     The  excess  energy  dissipated  in 
the  curve  is  obtained  from  the  backwater,  or  storage,  reach  upstream, 
where  there  is  a  compensating  reduction  in  the  total  energy  dissipated 
on  the  stream  bed.    At  some  distance  below  the  curve,  the  flow  will 
again  return  to  normal  conditions.     Thus,  for  straight  approach  and 
exit  reaches  of  sufficient  length,  there  would  be  no  net  excess  energy 
loss  through  the  curve.     This  is  a  general  rule  which  may  be  stated  for 
any  obstruction  or  resistance  in  an  open  channel. 

Therefore,  a  major  difficulty  in  determining  the  excess  rate  of 
energy  dissipation  due  to  a  curve  lies  in  the  definition  of  the  total 
reach  over  which  the  loss  is  measured.     In  order  to  make  open  channel 
tests  comparable  to  investigations  on  pipe  bends,  the  best  definition 
should  probably  be  based  on  a  reach  starting  above  the  curve  at  the 
point  of  maximum  specific  energy    and  extending  downstream  to  the  point 
at  which  uniform  flow  is  re-established.    The  exact  location  of  the 
downstream  point,  as  long  as  it  lies  in  the  region  of  uniform  flow,  is 
not  critical.    However,  since  the  upstream  point  of  measurement  is 
determined  by  the  backwater  curve,  it  is  clear  that  error  in  the  location 
of  this  point  will  modify  the  measured  results  accordingly. 

In  the  light  of  the  above  argument,  it  is  to  be  expected  that 
the  boundary  shear  stress  measured  at  Station  1  should  be  lower  than 
the  average  shear  stress,  yRSe,    computed  for  a  uniform  flow.     The  short 
length  of  the  test  channel  makes  it  impossible  to  evaluate  yRSe  accurately, 
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but,  as  shown  in  Table  II,  the  measured  values  of  t0  (at  Station  l)  are 
consistently  lower  than  the  estimated  average  shear  stress,  yRSe. 

It  is  seen  from  Figure  3U,  in  which  tx/ti  is  equivalent  to  the 
ratio,  S6x/Se]_,  that  the  energy  dissipation  is  only  moderately  increased 
in  the  upstream  portion  of  the  curve.    Near  the  curve  exit,  however,  the 
shear  stresses  rise  sharply  with  the  shift  of  momentum  to  the  outer  bank, 
and  these  high  shears  decrease  only  slowly  towards  the  normal  value  in 
the  channel  section  below  the  curve.     In  preparing  the  plots  to  determine 
U  for  the  tests  at  lower  curvature,  Runs  5>-7,  it  was  found  that  within 
the  early  part  of  the  curve  the  average  shear,  txAi,  decreases  to  a 
value  of  about  0.92,  before  increasing  towards  the  exit,    The  implication 
from  this  is  that,  for  decreasing  curvature,  the  storage  reach  may  extend 
further  into  the  bend.     The  resistance  of  the  curve  lies  primarily  in  the 
zone  in  which' the  irrotational  pattern  breaks  down.    Thus,  as  long  as  the 
flow  through  the  curve  takes  place  without  separation  and  is  essentially  in 
accordance  with  an  irrotational  pattern,  average  energy  losses  remain 
moderate  within  the  curve  itself.     The  greatest  part  of  the  loss  occurs 
below  the  point  of  breakdown,  over  a  reach  extending  well  below  the  curve. 

The  relative  amounts  of  the  loss,  which  occur  within  the  curve 
proper  and  over  the  downstream  reach,  vary  with  the  geometry  of  the  bend. 
Ito  (21)  has  shown  that,  in  pipe  bends  of  large  curvature,  as  much  as 
90°/6   of  the  bendhead  loss  may  take  place  in  the  downstream  tangent,, 
However,  for  moderate  curvatures  (d/rc  ^0.2),  only  about  I4.O  °/6   of  the 
energy  is  dissipated  in  the  sequent  reach,  due  to  the  reduction  in  the 
disturbance  generated  by  the  curving  boundary, 

The  general  trend  indicated  in  Figure  3£  suggests  an  increase  in 
excess  dissipation  with  increasing  curvature c    Because  of  the  small 
number  of  tests  available  for  this  comparison  and  the  relatively  large 
scatter  of  the  points,  no  attempt  was  made  to  give  concise  formulation  of 
this  trend.     Intuitively  such  a  trend  is  quite  plausible;  at  zero  curva- 
ture, the  resistance  is  that  of  a  straight  channel,  tl/t]_  =  1,  while  at 
maximum  curvature  the  resistance  should  attain  a  maximum  (finite)  value. 

Previous  investigations  on  losses  through  bends  in  pipes  (2),  as 
well  as  in  open  channels  (26),  have  revealed  that  there  may  occur  an 
anomaly  in  the  rate  increase  of  bend  resistance  with  increasing  curvature. 
For  bends  of  constant  roughness  and  deflection  angle  it  appears  that, 
with  increasing  curvature,  there  exists  an  intermediate  condition  of 
relative  minimum  resistance.     The  critical  value  of  curvature  at  which 
this  anomaly  occurs  has  not  been  clearly  defined,  nor  has  its  cause  been 
adequately  explained.     It  is  recognized  that  the  deviation  of  Runs  1  and 
2  from  the  general  trend  of  increasing  resistance  may  be  an  indication 
of  this  phenomenon  in  the  open  channel  tests.    However,  it  is  felt  that, 
because  of  the  small  amount  of  available  data,  there  is  no  justification 
in  this  case  for  assuming  such  a  discontinuous  behavior. 

In  concluding  this  discussion  it  can  be  stated  that  the  mean  loss 
coefficients  determined  for  flow  through  curves  from  the  usual  hydraulic 
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measurements  give  an  inadequate  indication  of  the  maximum  intensities  of 
local  shear  stresses.    For  flow  through  single  curves,  the  maximum  local 
shears  were  found  to  exceed  the  mean  initial  shears  by  over  100  while 
mean  sectional  shear  stresses  were  only  of  the  order  of  35>%  larger 
than  the  mean  stresses  at  the  entrance  of  the  curve.    For  conditions  of 
increasing  curvature ,  there  is  an  increasing  trend  in  the  mean  sectional 
shears,  as  well  as  in  the  total,  relative  resistance  of  the  entire  curved 
reach. 


E.     Scour  Patterns  in  Alluvial  Streams 

In  the  previous  section  the  energy  dissipated  by  a  curving  stream 
was  discussed  in  terms  of  an  average  resistance  over  the  total,  curved 
reach.    While  this  approach  is  of  interest  in  considering  the  role  of 
boundary  geometry  in  the  overall  energy  balance  of  the  stream,  it  tends 
to  obscure  the  local,  conditions  controlling  the  scour  pattern.    From  a 
comparison  of  the  shear  maps  with  the  mean  values  shown  in  Figures  3U 
and  3^5  it  is  clear  that,  in  an  erodible  channel  localized  areas  of  scour 
should  be  expected  at  various  sections  for  which  the  mean  shear  stress 
is  not  critically  high. 

For  the  understanding  of  scour  patterns,  therefore,  it  is  probably 
more  logical  to  consider  the  variation  in  magnitude  of  some  parameter  of 
local  shear  as  a  function  of  bend  geometry.    In  Figure  36,  the  magnitude 
of  the  maximum  measured  shear  stress  for  Runs  1-7  is  plotted  against 
relative  curvature,  w/rc.    Because  of  the  statistical  uncertainty  of 
actually  determining  an  extreme  value  in  these  tests,  the  quantitative 
validity  of  Figure  36  is  subject  to  question.    However,  the  trend  is 
clear,  and  it  may  be  concluded  that,  in  an  erodible  trapezoidal  channel, 
these  would  represent  points  at  which  erosion  would  first  appear.  As 
noted  in  the  earlier  discussions,  with  increasing  curvature  this  point 
of  initial  scour  moves  to  a  position  further  upstream.    On  the  basis  of 
these  tests,  this  shift  is  apparently  quite  abrupt,  occurring  at  a 
curvature  in  the  range  O.67  <  w/rc  <  0.73. 

An  additional  parameter,  which  may  also  be  relevant  in  an  approach 
to  the  problem  of  channel  stability  and  the  design  of  revetments,  is  the 
variation  in  areal  extent  of  a  zone  subjected  to  shears  greater  than  some 
critical  magnitude.    Consider,  for  example,  the  zones  in  the  curved  reach 
(Runs  1-7)  enclosed  by  the  contour  t0/       =  l.U.     These  areas,  expressed 
as  a  fraction  of  the  area  of  the  test  reach,  are  plotted  in  Figure  37* 
As  might  be  expected,  the  relative  area  exposed  to  shears  equal,  to  or 
greater  than  a  given  value  increases  with  increasing  curvature.  For 
channels  of  fixed  boundary  geometry  (Runs  1-h,  or  5>-7),  the  increase  in 
depth  and  curvature  at  flood  conditions  will  lead  to  scour  which  is 
increased  both  in  extent  and  intensity. 

Before  it  becomes  feasible  to  apply  these  experimental,  results 
directly  to  natural  streams,  further  information  is  needed  on  the  effects 
of  boundary  deformation  on  the  pattern  of  scour.     It  is  evident  that,  at 
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a  point  of  high  local  shear  stress ,  scour  will  occur.    It  cannot  be 
predicted,  however,  to  what  degree  the  resulting  boundary  depression  will 
modify  the  shear  either  at  the  point  or  in  the  surrounding  area.  In 
addition,  the  deposition  which  may  be  expected  in  the  zones  of  reduced 
shear  should  further  modify  the  flow  pattern  by  localized  aggradation  of 
the  boundary. 

The  variation  in  erodibility  of  sorted  sediments  adds  still 
another  complicating  factor  to  this  problem.    Sundborg  (U3)  nas  shown 
that  the  sorting  of  a  heterogeneous  bottom  sediment  depends  on  the  size 
range  as  well  as  on  the  particle  shapes  and  the  nature  of  the  sediment 
mixture.     In  some  sediments,  the  fine  material  is  winnowed  out,  leaving 
an  erosion  pavement  which  can  withstand  much  higher  shears.    Under  other 
conditions  the  residual  sediment  is  comprised  of  the  finer  material, 
which  may  then  develop  into  ripples  and  thus  alter  the  effective  bottom 
roughness.    Valuable  information  could  be  gained  from  both  laboratory  and 
field  work  in  which  systematic  studies  were  made  to  determine  the  sorting 
of  bottom  sediments  under  varied  shear  stresses.    Of  particular  interest 
to  the  present  study  would  be  data  on  the  size  distribution  across 
transverse  sections  in  stream  bends.    As  far  as  could  be  determined,  this 
work  has  never  been  done. 

In  spite  of  these  .still  unsolved  problems,  a  quite  logical  rela- 
tionship can  be  observed  between  the  scour  patterns  in  alluvial  streams 
and  the  boundary  shear  distributions  measured  in  the  channel,  experiments. 
Figure  38~x~  shows  a  short  reach  of  Buffalo  Creek,  Erie  County,  New  York, 
in  an  area  where  the  stream  flows  through  an  unconsolidated  mixture  of 
glacial  sands  and  silty  clays.    The  stream  discharge  of  2, £00  cfs  to 
which  the  map  corresponds  is  based  on  the  approximate  conditions  of  annual 
spring  flooding. 

By  comparison  of  the  section  profiles  with  the  map,  the  scour 
pattern  through  the  reach  can  be  observed.    Before  each  bend  the  stream 
assumes  an  approximately  trapezoidal  shape,  indicating  a  symmetrical 
flow  distribution.    The  tendency  for  the  flow  to  cross  to  the  outer  bank 
under  the  influence  of  the  helicoidal  motion  extends  quite  far  upstream 
in  each  bend,  as  evidenced  by  the  early  dislocation  of  the  thalweg. 
Between  the  second  and  third  bends,  an  oblique  bar  is  evident  at  the 
crossing  point,  with  a  discontinuity  of  the  thalweg.    This  is  a  type  of 
bar  configuration  which  is  common  to  many  sinuous  streams.    The  scour 
extending  down  to  the  inner  bank  of  Bend  3  is  due  to  the  tendency  to 
irrotational  motion  of  the  flow  entering  the  curve,  as  well  as  to  the 
asymmetry  of  the  flow  leaving  Bend  2.    A  similar  condition  may  be 
observed  for  the  reverse  curve  tests  in  the  channel,  Figures  30*  31»  ^n 
these  bends  of  moderate  curvature,  the  helicoidal  motion  should  be 
developed  early  in  the  ^curve,  and  the  deeper  channel  on  the  concave  bank 
is  thus  extended  quite  far  upstream. 

 In  closing  the  discussions  of  these  experiments  it  may  be  remarked 

"This  material  was  supplied  by  Messrs.  D.  A.  Parsons  and  R.  P.  Apmann,  of 
the  Agricultural  Research  Service. 


-66- 


that  the  magnitudes  of  the  measured  shears,  which  might  be  expected  to 
cause  bank  instability,  compare  favorably  to  field  observations  on  erosion 
in  curves.    Lane  (22)  set  forth  stability  criteria  based  on  tractive 
force  studies  conducted  in  unlined  channels.    For  canals  ranging  in 
curvature  from  straight  to  very  sinuous,  the  stability  criterion  is  given 
as  a  fractional  reduction  in  the  permissible  tractive  force.    It  is  stated 
that,  for  very  sinuous  canals,  the  permissible  tractive  force  (equivalent 
to  tq)  is  reduced  by  a  factor  of  0.6.     Thus,  the  maximum  relative  tractive 
force  to  be  expected  in  such  curves  is  t0/       =  1/0.6  =  1.7. 

The  degree  of  curvature  as  defined  by  "very  sinuous"  corresponds 
to  the  areal  configuration  "typical  of  canals  in  foothills  or  mountainous 
topography"  [(22),  p. 12^7],    A  reasonable  estimate  of  curvatures  for  such 
streams  appears  to  be  of  the  order  of  0.U  <  w/rc<06.    While  the  maximum 
shears  measured  in  the  channel  are  somewhat  greater  than  this,  the  agree- 
ment is  quite  good. 

It  may  be  concluded  that  the  local  tractive  forces  measured  under 
these  simplified  channel  conditions  do  provide  a  useful  indication  of  the 
distribution  and  intensity  of  erosive  attack  on  the  boundaries  of  alluvial 
channels.    In  order  to  establish  closer  correlation  between  the  laboratory 
and  field  conditions,  further  effort  should  be  directed  towards  two  major 
areas:     (l)  observation  techniques  leading  to  more  concise  quantification 
of  field  data  and    (2)  the  interaction  between  flowing  water  and  the 
stream  sediments.    The  latter  problem  includes  the  role  of  entrained 
sediment  as  it  modifies  the  dynamical  properties  of  the  transporting 
fluid,  as  well,  as  the  effect  of  boundary  alteration  arising  from  movement 
of  the  bed  material. 
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VI.  CONCLUSIONS 


The  results  of  this  investigation  of  flow  through  curved  trape- 
zoidal channels  permit  a  number  of  important  conclusions. 

A.    Flow  Through  Single  Curves 

1.  For  channels  composed  of  a  single  60-degree  curve  with 
straight  approach  and  exit  tangents ,  the  flow  pattern  within  the  curve 
tends  to  conform  to  that  of  a  free  vortex.    Boundary  layer  effects, 
combined  with  the  centripetal  acceleration  in  the  curving  flow,  cause 
a  trend  to  separation  along  the  concave  bank  near  the  curve  entrance 
and  a  more  pronounced  separation  zone  originating  at  the  convex  bank 
towards  the  exit  of  the  curve.     In  addition  to  separation,  the  fric- 
tional  effects  produce  the  expected  helicoidal  motion. 

2.  The  two  frictional  aspects  of  the  flow,  separation  and  the 
secondary,  helicoidal  current,  tend  to  produce  a  breakdown  in  the 
irrotational  type  of  motion,  with  the  ultimate  establishment  of  high 
velocities  along  the  concave  bank.     The  section  of  the  curve  at  which  the 
stream  momentum  is  transferred  from  the  inner  to  the  outer  bank  varies 
with  the  bend  geometry.    This  crossover  point  occurs  earlier  in  the  curve 
for  streams  of  decreasing  curvature.    Under  conditions  of  greatest  curva- 
ture, the  irrotational  type  of  pattern  prevails  throughout  the  curve 
itself f  and  at  the  exit  section  the  crossover  is  abrupt,  due  to  the  large 
separation  zone  on  the  inner  bank. 

3.  The  distribution  and  the  magnitudes  of  the  boundary  shear 
stresses  produced  by  the  curved  flow  cannot  at  present  be  predicted 
quantitatively,  either  by  theory  or  by  correlation  with  measured  velo- 
city distributions.     In  general,  the  locations  of  the  shear  maxima 
were  found  to  be  associated  with  the  course  of  the  filament  of  highest 
velocity  and  with  the  zones  of  locally  accelerated  motion.    Thus,  high 
stresses  are  found  along  the  inner  part  of  the  curve    and  near  the  outer 
bank  below  the  curve  exit,  with  a  strongly  asymmetrical  shear  distri- 
bution persisting  for  a  considerable  distance  in  the  downstream  tangent 
reach. 

U.     The  location  and  the  relative  intensity  of  the  maximum  local 
shear  stress  in  a  curve  are  controlled  primarily  by  the  stream  geometry. 
For  streams  of  increasing  curvature,  the  relative  magnitude  of  the  shear 
maximum  increases,  and  it  tends  to  occur  earlier  in  the  curved  reach. 
For  the  range  of  curvatures  studied  in  the  smooth  channels,  0.29  <  w/rc  < 
0.8,  the  shear  maximum  increased  from  t0/  t©  =  1.6  to  2.U. 

5>.     The  total  head  loss  for  the  curve,  as  well  as  the  average 
rates  of  loss  at  sections  throughout  the  curved  reach,  may  be  computed 
from  the  measured  shear  distributions.    This  one-dimensional  treatment 
of  the  flow  gives  inadequate  indication  of  either  the  intensities  or  the 
peripheral  locations  of  the  highest  shears;  however,  it  demonstrates 
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that  the  excess  stream  energy  dissipated  due  to  a  curve  increases  with 
increasing  curvature. 

6.    The  introduction  of  boundary  roughness  does  not  appear  to 
modify  the  shear  pattern  within  the  curve  itself.    However ,  in  the 
downstream  reach  below  the  crossover  point,,  the  shears  are  substantially 
higher ,  due  to  the  greater  development  of  the  helicoidal  motion  over 
the  rough  boundary. 

B.     Compound  Curve  Systems 

1.  The  technique  of  simulating  the  effects  of  upstream  curva- 
ture by  artificially  reproducing  the  related  velocity  distributions  in 
the  approach  channel  has  proved  to  be  useful  and  economical  for  estab- 
lishing shear  patterns  for  compound  curves.  The  quantitative  validity 
of  this  method  remains  subject  to  detailed  confirmation. 

2.  The  effect  of  a  reversed  curvature  is  simulated  with  high 
velocities  near  the  inside  bank  at  the  entrance  section.     This  condition 
reinforces  the  free  vortex  pattern ,  increasing  the  superelevation  to  a 
slight  degree.    High  shear  stresses  are  produced  along  the  inner  bank. 
Separation  towards  the  exit  section  is  prevented,  and  thus  a  more  sym- 
metrical flow  condition  exists  in  the  downstream  tangent  with  shear 
stresses  much  reduced,  as  compared  to  the  case  of  the  single  curve,, 

3.  With  high  velocities  produced  near  the  outside  bank  at  the 
entrance  section  a  double  curve  configuration  is  simulated  with  a  short 
straight  section  between  the  curves.     This  approach  condition  suppresses 
the  free  vortex  pattern  in  the  curve;  a  larger  separation  zone  is  formed 
near  the  inside  bank;  and  higher  shear  stresses  are  encountered  along 
the  outer  bank  throughout  the  test  reach  than  for  the  single  curve. 

U.    In  general,  it  may  be  concluded  that  longer  curves  will  not 
materially  increase  maximum  shear  intensities;  by  contrast,  reversal  of 
curvature  in  compound  systems  leads  to  a  substantial  increase  in  local 
shear  intensity.     The  absolute  maximum  relative  shear  stress  was  observed 
to  be  3-0  in  the  simulated  reverse  curve. 


C.     General.  Observations 

1.  For  stream  bends  of  moderate  curva  ture,  the  superelevation 
of  the  water  surface  is  quite  insensitive  to  variation  in  the  radial 
distribution  of  velocities    and  may  be  predicted  adequately  by  a  one- 
dimensional,  treatment  of  the  flow.     The  total  superelevation  appears  to 
be  a  function  only  of  the  mean  momentum  of  the  stream    and  the  relative 
curvature  of  the  bend. 

2.  For  a  given  channel  alignment,  the  relative  shear  patterns  in 
a  curve  are  in  general,  not  greatly  modified  by  variations  in  depth  and 


-68- 


velocity  at  the  curve  entrance,  but  seem  to  depend  primarily  on  the 
boundary  geometry.    For  the  range  of  flow  conditions  which  can  occur 
in  a  trapezoidal  channel,  the  intensities  and  the  areal  extent  of  the 
highest  shears  increase  with  increasing  depth  and  curvature. 

3.     This  study  has  demonstrated  the  importance  of  boundary 
shear  stress  patterns  with  respect  to  the  ultimate  understanding  of 
erosion  and  deposition  processes  in  channel  curves.    It  has  been  shown 
that  the  trends  in  the  measured  shear  distributions  conform  quite  closely 
to  the  scour  patterns  observed  in  sinuous  alluvial  streams,  and  that, 
insofar  as  the  limited  field  data  permit  a  comparison,  the  magnitudes  of 
the  measured  stresses  suggest  agreement  with  the  intensity  of  the  channel 
scour. 

U.    While  the  orientation  of  the  shear  stresses  was  neglected 
in  this  study  and  is  not  believed  to  affect  the  stated  results,  it  may 
be  concluded  that  it  assumes  considerable  importance  with  respect  to 
bank  slope  stability  in  channels  through  uniform  non-cohesive  materials. 

5.    The  distribution  of  the  shear  stresses  over  the  boundary  is 
conveniently  determined  from  surface  Pi tot  tube  measurements .  This 
technique  as  previously  developed  by  Preston  for  air  flows  has  been 
successfully  adapted  to  the  present  problem  for  smooth  as  well  as  rough 
boundary  surfaces. 
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VIII.  APPENDIX 


Superelevation  in  Curved  Flow 


Consider  the  flow  through  a  circular  curve  in  a  rectangular 
channel  as  shown  in  the  accompanying  sketch.    Neglecting  entrance 
transition  effects ,  let  it  be  assumed  that  all  motion  is  concentric 


and  that  there  is  no  frictional  dissipation  of  energy  on  the  stream 
boundaries.    Thus,  the  mean  specific  energy  is  constant  at  all  sections 5 
and  may  be  computed  from  the  approach  conditions: 

r  2 
H  =   - — -      f    Hdr-H    =y    +  ~-  .  (A-l) 

r. 
i 


Since  the  motion  of  each  fluid  particle  is  circular }  equilibrium  of 
radial  forces  is  given  by: 

dp_  _  pu(r)2  / 
dr  r 


and,  for  the  free  surf ace,  the  superelevation  becomes: 


However,  the  computed  flow  through  a  curve  in  a  stream  of  specified  dis- 
charge must  be  consistent  with  the  requirements  of  the  conservation  of 
mass  and  energy.     Thus  the  assigned  function  u(r)  and  the  depth  within 
the  curve  y(r)  must  satisfy: 


H  =  H  = 
o 


r   -  r . 

o  1 


r 

f\0 


Wr)t2^)dr, 


(A-U) 


and 


Q  =  w  yQ  V  = 


y(r)  u(r)  dr 


(A-5) 


As  a  first  model  for  curved  flow,  the  superelevation  is  computed  for  the 
condition  of  specific  energy  constant  in  the  radial  direction.  Then: 

H(r)  =  y(r)  +  =  Hq  ,  (A-6) 

and 

¥=~~lr  ■  (A-7) 

dr      g  dr 

From  equations  (A-2)  and  (A-7),  the  velocity  distribution  function  may 
be  computed  as  that  for  a  free  vortex: 

u  =  J:  ,  (A-8) 

and  the  total  superelevation  is  given  by: 

K  2 

2g  r .  2  r 
1  o 
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As  shown  in  the  definition  sketch  the  effective  mean  radius , 
r  9  is  defined  such  that: 


u 


(r  =  rj  =  7' 


Hence, 


K-  =  V  r, 
J.  " 


(A-10) 


and 


Vr 


u  - 


Then  the  mean  depth  within  the  curve  may  be  obtained  as: 

r 

o  f\0 

,2 


^m  r 


r. 
1 


y(r)dr  a 


r  -  r. 
o     i  \J 
r. 
l 


(HQ-^_)dr,  (A-ll) 


where  y(r)  is  as  defined  implicitly  in  equation  (A-6) .  The  resulting 
expression: 


V 


2  r, 


y   =  H   —  , 

m       o    2g  r  r.  ' 


may  be  arranged  as: 


m 


(r„/r  )2 
»  c 

1-  (— ): 

c 


(A-12) 


Equation  (A-12)  indicates  that  the  mean  depth  in  the  curve  may  not  be 
assumed  a  priori  eaual  to  the  depth  in  the  approach  flow. 


Equation  (A-£)  becomes: 
r 


Q 


-f 


V3r  3  r.2  -r  2 


r. 
i 


— (H  --5— ^)dr  =  H  Vrv,  ln(-g)  +-1JL(  1  ,    °J  ,(A-13) 
r     o    2g  r2'  o     #       r.        ag      r  2  r.2'  5V 

&  1         6  01 


which  leads  to  the  cubic  equation: 
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r  r 


0         In'  c 


c 


1- 


w 
2r 


7 


2[1-  (-^-)2] 

 c 

F2 


=  0  (A-1U) 


Then,  equations  (A-9) ,  (A-10) ,  and  (A-lU)  may  be  combined  to  give: 


y2  y2 

(Mj  .(SSL) 


(2  +P2)  In 


 c 

&-<£■)  ] 

c 


=    0  (A-15) 


Equation  (A-l£)  is  plotted  in  Figure  A-l  showing  the  relative 
superelevation,  hy/y.o,  as  a  function  of  Froude  number  for  three  values 
of  relative  curvature,  w/rc.    As  shown  in  the  text  of  this  report  (see 
pp.  U8  ff.)  the  initial  assumption  of  constant  mean  depth  in  the  curve 
leads  to  the  expression: 


2g  r 


1-  (— ) 


(lib) 


or 


o  c 


(A-16) 


which  is  also  plotted  in  Figure  A-l.    It  is  seen  from  the  figure  that  at 
higher  Froude  numbers  equation  (A-16)  gives  values  of  Ay/y0  which  are  too 
low, although  at  moderate  curvatures  the  differences  are  not  great  and 
the  assumption  ym  =  y0  is  reasonable. 

As  the  second  flow  condition,  it  is  assumed  that  the  velocities 
increase  linearly  in  the  radial  direction.    This  flow  distribution,  which 
is  seen  near  the  end  of  a  long  river  bend  (cf  Runs  £,6,7,10 ,11) corres- 
ponds to  a  forced  vortex.    Under  the  assumption, 


u 


K2  r  , 


(A-17) 
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the  total  superelevation  is  given  by: 


K 

2g     o  1 


(A-18) 


Again  introducing  the  effective  mean  radius ,  r„,  the  expression  reduces 
to: 


r  r 
c 


(A-19) 


By  the  assumption  of  constant  mean  specific  energy,  the  mean  depth 
in  the  curve  is  given  by: 


V    1    ro3  -  ri3 


ym       o  "  2g  r„s  3(r  -  r  ) 

A  O  1 


(A-20) 


which  may  be  expressed  by: 


ym  _ 

^  =  1  ~ 


(-) 


3 


(A-21) 


Now  in  order  to  evaluate  the  continuity  relation,  (A-£),  it  is 
necessary  to  further  define  y(r),  since  in  this  case  the  function,  H(r), 
is  an  unknown.    Expressing  y(r)  in  terms  of  the  depth  at  the  inner  bank: 


y(r)  =  7T.  +  6  y(r)  =  yr 
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2,  r  dr, 


(A-22) 


and 
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no 


r   -  r . 

o      1  u 


r 

AO 


y(r)dr  = 


Jri    r   - r 


o      i  0 


r 
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i    r  • 


v* 


r  2ct 


r  dr 


dr 
(A- 23) 


This  equation  can  be  integrated  and  rearranged  to  give: 
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y  r 

...  =  IS  .  P2(_£) 

y  r 
o      J  o  * 
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■(3-—) 
2rcU  2rc; 


(A-2U) 


-79- 


Introducing  equation  (A-22)  into  the  continuity  expression,  equation 
(A-£)  becomes: 


r 

no 
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r . 
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*\,2g 


r  dr 
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r  dr  . 


(A-25) 


Equations  (A-19),  (A- 21) ,   (A-2U),  and  (A-2£)  may  then  be  combined  to 
give : 

-#)(T-)(2+  r2)  +  m/2   F(^-)     =0  (A-26) 


3/2 


3-  (— ) 
V2rc; 


y  r 
17  o  c 


Equation  (A-26)  is  plotted  in  Figure  A-2,  together  with  the  expression 
based  on  the  assumption  of  constant  mean  depth: 

V2  2w 


which  may  be  rewritten  as: 


It  may  be  seen  from  the  figure  that,  except  for  conditions  of  severe 
curvature,  the  two  expressions  give  virtually  identical  values  of  Ay/y0, 
and  that  the  assumption  of  constant  mean  depth  is  completely  justified. 

From  the  foregoing  developments  which  lead  to  equations  (A-l£) 
and  (A-26)  it  will  be  seen  that,  in  general,  if  the  velocity  distribution 
for  flow  in  an  open  circular  curve  can  be  expressed  as  some  integrable 
function: 

u  =  u(r)  , 

then  an  exact  solution  for  the  total  superelevation  may  be  found,  which 
is  based  on  the  single  fundamental  assumption  of  constant  total  energy. 
Although  the  complete  form  of  the  transverse  surface  profile  has  not  been 
given  here,  it  can  be  shown  that  the  explicit  function  for  the  depth  at 
any  radius,  equation  (A-22),  can  be  obtained  through  solution  of  expressions 
of  the  form  of  equations  ( A- 21) , (A-23) ,  and  (A-2U). 

It  is  again  recalled  that,  consistent  with  the  above  assumption, 
no  restrictions  need  be  placed  on  the  geometry  of  the  water  surface.  For 


1  + 


w^ 
12r  2 

r>  —i 


(ll-o) 


(A-27) 
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example,  equations  (A-12)  and  (A-21)  indicate  that  an  a  priori  assumption 
of  ym/y0  =  1  is  not  justified.    However ,  for  a  broad  range  of  curvatures 
and  Froude  numbers,  essentially  covering  normal  stream  conditions, 
Figures  A-l  and  A- 2  indicate  that  the  error  introduced  by  this  simpli- 
fying assumption  is  in  fact  negligible. 
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